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Abstract 


The development of new multidimensional NMR techniques and methodologies, in recent 
years, has facilitated in obtaining high precision protein structures in solution. However in the 
case of large size proteins (M r > 20 kDa), poor resolution due of severe spectral overlap 
comphcates their analysis. In order to address this problem, we have proposed a novel 
methodology called amino acid selective "unlabeling" . In this approach, one or more specific 
amino acid residues in a protein are selectively "unlabeled” while uniformly or fractionally Li C 
or/and N labeling the rest. The cross peaks belonging to the unlabeled amino acid residues are 
thus absent in various NMR spectra recorded with such a protein sample, leading to significant 
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assignments and stereospec.fic resonance assignments of diastereotopic methyl groups of Val 
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shifts (PCS) and residual dipolar coupling (RDC) in a 15 kDa calcium binding protein from 
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Introduction 

The development of new experimental techniques and methodologies in the field 
of nuclear magnetic resonance (NMR) spectroscopy, in recent years, has facilitated in 
determining the 3D structures of proteins in solution to high precision 1 ' 3 . The size of 
proteins that are currently being analyzed by NMR is increasing steadily and attempts 
are on to elucidate structures of proteins as large as 1 1 0 kDa 4 . However, application of 
this technique has been limited, in large molecular weight proteins, by the relaxation 
properties of various nuclei involved in any experiment and the resolution of spectra. 
In many cases, severe spectra overlaps hamper analyses even while using 
multidimensional experiments combined with 13 C/ I5 N/ 2 H isotopes enrichment of the 
protein under investigation. Thus, one of the present challenges in protein NMR today 
is to identity and provide solutions to reduce overlaps in NMR spectra of large 
proteins. 

A novel methodology, which addresses this problem of spectral overlap, is the so 
called amino acid selective "unlabeling" 5,6 . This is based on selective unlabeling of 
amino acid residues in a uniformly or fractionally 13 C or/and 15 N labeled protein, 
which simplify the multi-dimensional heteronuclear NMR spectra. A particular amino 
acid in a protein is 'unlabeled' by feeding the host microorganism expressing the 
protein of interest, with 15 N labeled ammonium chloride or/and I3 C labeled glucose as 
the sole source of nitrogen and carbon, respectively, along with the desired amino 
acids(s) to be unlabeled, in unlabeled form. This renders the desired amino acid 
residue(s)in the protein unlabeled, as a result of which cross peak due to these residues 
are not observed in any of the double and triple resonance spectra. A comparison of 
such a spectrum with that of a control spectrum involving uniformly l5 N or/and l3 C 
labeled protein enables one to distinguish peaks, and thereby the corresponding 
chemical shift of nuclei belonging to the unlabeled amino acid residues(s). 

The utility of this methodology has been demonstrated previously for residue 
specific assignments 5 and stereospecific assignments of diasterotopic methyl groups in 
Val and Leu residues in large proteins 6 . In the former case, amino acid selective 
unlabeling combined with an automated approach, TATAPRO 7 , was shown to aid in 
sequence specific resonance assignments, which form the first step towards unraveling 
3D structures of proteins by NMR 1 . On the other hand, the second application 
involved selective unlabeling of amino acid residues in a fractionally (10%) Re- 
labeled protein, to obtain unambiguous stereospecific assignments of diastereotopic 
methyl groups in Val and Leu residues 6 . Such stereospecific assignments provide 
crucial NMR constraints in computing and refining the 3D structure of proteins 8 . 
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In this paper, we demonstrate yet another utility of this method, namely, the 
measurement of pseudo-contact shifts (PCS) and residual dipolar couplings (RDC) for 
proteins with paramagnetic centers. PCS refers to shifts in resonance frequencies 
(chemical shifts) of nuclei arising due to their hyperfine interaction with an unpaired 
electron residing on the paramagnetic center. In a metal centered approximation, PCS 
is given by 9 : 


8 pc (pp m ) = 


1 

12^r 3 





( 1 ) 


where AXax and Ax* are the axial and rhombic magnetic susceptibility anisotropy 
induced by the paramagnetic center and /„ m, and n, refer to the direction cosines of 
the position vector of atom i (r,) with respect to the orthogonal system formed by the 
principle axes of the magnetic susceptibility tensor y? ara . As evident from this 
equation, PCS contains structural information on both the distances and angles in the 
frame of the magnetic susceptibility tensor. Moreover, the distance dependence is of 
the 1/r 3 type, which, when compared to relaxation rate enhancements (1/r 6 ) type, 
decays more slowly with distance and thus can provide information at a greater 
distance from the metal ion than those obtained via the conventional NOESY 
experiments. 

In addition, the paramagnetic centers, which are characterized by high anisotropy 
in magnetic susceptibility, induce a partial orientation of the molecules in high 
magnetic field (B 0 ). This is owing to the fact that different orientations of the protein 
(with a paramagnetic center) with respect to B 0 give rise to different interaction 
energies, resulting in a residual alignment with respect to B 0 . Thus, the dipolar 
coupling between two nuclear spins do not average out to zero, as it usually happens 
in liquids, giving rise to residual dipolar coupling (RDC). These RDCs contain 
information on the orientation of the vector connecting the two dipoles with respect to 
the molecular magnetic susceptibility tensor (x“ o1 ) as given by 9 : 

Ax**- (3cos 2 0 - 1)+ |Ax M (sin 2 6cos2<|») 


i Bq y hK n ^ 2 
An 15 kT 4rc 2 rl N 


RDC(Hz) = 


( 2 ) 



276 


H. L. ATREYA et al. 


The strength of such individual RDCs gives an estimate of the degree of alignment, 
which can then be used as a form of structural constraint in restrained molecular 
dynamics calculations, for protein structure determination 9 . 

In the present study, involving a calcium binding protein from Entamoeba 
histolytica (EhC aBP), lanthanum ions were chosen as paramagnetic probes. This is 
owing to the fact that lanthanide ions have similar ionic radii as calcium and thus can 
replace the latter in a protein. Further, they posses large anisotropic magnetic 
susceptibilities, a property very useful for obtaining long range structural information 
as discussed above. However, in the case of lanthanum loaded E/zCaBP, severe 
spectral overlaps and large resonance line widths prevent accurate measurement of 
both PCSs, and RDCs. This promoted us to adopt the method of selective amino acid 
"unlabeling", which resulted in significant spectral simplification and facilitated in 
accurate measurement of PCSs and RDCs. The following paragraphs describe the 
outcome of such an endeavor. 


Materials and Methods 


(i) Sample preparation : 

Recombinant proteins are enriched isotopically with lj C and 15 N by growing the 
microbial host in a M9 salt medium, supplemented with a l5 N labeled ammonium 
chloride as the sole source of nitrogen and i3 C labeled glucose as the sole source of 
carbon 10 . £/zCaBP was cloned and over-expressed in E. coli BL21 (DE3) strain 
containing pET-3 c expression system. Uniformly 15 N labeled EhCaBP and amino acid 
selectively unlabeled-' ^ labeled EhCaB? samples were obtained using the protocol 
described earlier 5,7 . In the present study, two samples were prepared : (i) uniformly 
i5 N labeled ETzCaBP (hereafter referred to as Sample 1) and (ii) Arg, Phe, lie, Leu and 
Lys unlabeled-uniformly 15 N labeled EhC aBP (hereafter referred to as Sample 2). The 
choice of the five amino acids was based on the fact that, these residues are abundant 
in the protein and at the same time away from the metal binding sites. 

(ii) NMR experiments : 

NMR experiments were carried out on a Varian Unity + 600 NMR spectrometer 
equipped with a pulsed-field-gradient unit and triple resonance probe with actively 
shielded Z-gradients, operating at ’H frequency of 600.051 MHz. All measurements 
were performed with 0.6 ml samples of approximately 2 mM EhCaS? in 30 mM 
CaCf and 50 mM deuterated TRIS buffer (pH = 6.0) and at a temperature of 35 °C, in 
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a mixed solvent of 90% H 2 O and 10% 2 H20. For the measurement of PCS, a 2D 
sensitivity enhanced [ 15 N-'H]-HSQC u was recorded on lanthanide substituted 
EhC aBP, while RDC measurements were obtained via co i ( 15 N- 1 H) coupled 2D [ 15 N- 
'H]-HSQC. Data transformation and processing were done on Silicon Graphics 
workstation (R10000 based Indio II Solid Impact Graphics) using the FELIX97 
software (MSI, U.S.A.) 

Dy 3+ was chosen as the lanthanide ion in the present study. A stock solution 
containing 20 mM DyCl 3 (obtained from SIGMA) was prepared. Addition of the salt 
to EhC aBP was carried out in the form of titration, wherein 5 pi of the stock solution 
(corresponding to 0.1:1 metal-protein ratio) was added at every step. PCS and RDCs 
were measured at 1:1, 2:1, 3:1 and 4:1 metal protein ratio. 

Results and Discussion 

EhCaBP is a 14.7 kDa monomeric polypeptide (134 amino acid residues) 
possessing four calcium- binding loops characteristic of EF-hand family of calcium 
binding proteins. This is evident from its primary sequence shown in Scheme I, where 
highly homologous Ca 2+ binding loop segments are highlighted : 

ME ALFKEID VNGDGA VSYEEVKAFV SKKRAIKNEQ 

LLQLIFKSIZM2)GyVG£'/DgA®FAKFYGSIQGQDLSDDKI 

GLKVL YKLMD VDGDGKL TKEEVTSFFKKHGIE 

KV AEQVMKADANGDGYITLEEFLEF S L 

Scheme I 

The percentages of various amino acids residues present in £/zCaBP are Ala- 
7.5%, Arg-0.7%, Asn-3.7%, Asp-9.7%, Cys-0%, Gln-4.5%, Glu-9.7%, Gly-9.0%, 
His-0.7%, Ile-6.7%, Leu-9.7%, Lys-11.9%, Met-2.2%, Phe-6.7%%, Pro-0%, Ser- 
5.2%, Thr-2.2%, Trp- 0%, Tyr-3.0%, Val-6.7%. Fig. la shows the sensitivity- 
enhanced 2D [^N-'^HSQC spectrum of uniformly 15 N labeled EhCaBP in its Ca 2+ 
bound form (Sample 1). Titration of [u- 15 N] labeled (Ca 2+ ) 4 - EhCaBP (holo) with Dy 3+ 
resulted in several changes in the HSQC spectrum due to the binding of Dy 3+ ions to 
EhCaBP. Several pseudo^contact shifted peaks were seen in the FISQC spectrum, with 
some of the cross peaks showing severe line broadening. Such an observation is 
illustrated in Fig. lb, which shows a sensitivity enhanced 2D [ 15 N- X H]HSQC spectrum 
of uniformly 15 N labeled EhCaBP at a metal (Dy 3+ ') : protein ratio of 1 : 1. As is 
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evident from this Fig., it is impossible to assign pseudo-contact shifted peaks, 
particularly those having ! H chemical shifts between 7.00 ppm and 9.5 ppm, 
precluding accurate measurements of PCSs and RDCs in this protein. 

In principle, it is possible to reduce such spectral overlap by recording 3D NMR 
experiments, involving one or more heteronuclei. However, it is important to note that 
in proteins containing paramagnetic centers, resonance lines are considerably 
broadened owing to strong nuclear-electron dipolar interactions 9 . Thus, due to the fast 
transverse relaxation of heteronuclei involved, 3D NMR experiments are considerably 
reduced in sensitivity as compared to 2D experiments. One is then forced to use 2D 
[ l3 C/ 15 N-'H] HSQC spectrum only, to obtain measurements of PCSs and RDCs. In 
order to simplify this spectrum, we prepared a sample wherein the five amino acid 
residues, namely Arg, lie, Leu, Lys and Phe, were selectively unlabeled (Sample 2). 
The choice of these amino acids was based on two facts. First, these residues are 
among the most abundant in the protein and do not cross-matabolize to other amino 
acids during their biosynthesis 10 . Secondly, as evident from Scheme I, these amino 
acid residues occur less frequently in the four Ca 2+ -binding loops of EhC aBP as 
compared to other residues implying that, most of these amino acid residues are away 
from the Ca 2+ binding loops. Although, information on these amino acids are lost due 
to their selective unlabeling, we can measure PCSs and RDCs for the rest of the amino 
acid residues in the protein that lie in the vicinity of the metal binding sites. 

Selective unlabeling of these five residues resulted in the absence of cross peaks 
belonging to a total of 47 residues. Fig. lc and Id show the 2D [ I5 N-’H] HSQC 
spectrum recorded on Sample 2 in the Ca 2+ -bound form and at a metal (Dy 3+ ) : protein 
ratio of 1:1, respectively. A dramatic improvement in the resolution, obtained as a 
result of selective unlabeling, is evident from these figures, as compared to that of 
Sample 1 (Fig. la and lb). Cross peaks belonging to several residues in the vicinity of 
the metal binding sites are now amenable for accurate measurement of PCSs and 
RDCs as discussed in detail below. 

Fig. 2 A-C show selected regions of an uncoupled and coupled 2D [ l5 N-'H] 
HSQC spectrum recorded on Sample 1 and 2 at a Dy 3+ -protein ratio of 1 : 1 . Each of 
the Fig. 2 A-C is divided into five panels (A;, B, and C;, where i = 1-5). Panels 1 and 
2, show a selected region of spectrum extracted from an uncoupled 2D [ 15 N-’H] 
HSQC spectrum of Sample 1, while panels 3 and 4 show the same region for sample 
2. Panel 5, extracted from a 2D coupled [^N-'H] HSQC spectrum corresponds, in 
each case, to the region shown in panel 4 and is used in measuring RDCs for pseudo- 
contact shifted cross peaks. These regions serve as illustrative examples of 
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unambiguous measurements of PCSs and RDCs that are facilitated as a result of 
selective unlabeling. All the assignments shown in these figures are based on 
sequence specific resonance assignments of *H and 15 N spins in EhCaB? that have 
been reported earlier 12 . 



Fig. 1-2D sensitivity enhanced [ ls N-'H]HSQC spectrum of (a) uniformly 15 N-labeled (Ca 2+ ) 4 . 
EhCaB?, (b) uniformly ,s N-labeled EhCaB? at Dy 3+ : protein ratio of 1 : 1, (c) Arg, 
lie, Leu, Lys and Phe unlabeled-rest unlabeled 15 N labeled (Ca 2+ )„. EhCaB? 
and (d) Arg, lie, Leu, Lys and Phe unlabeled-rest 1S N 
labeled EhCaB? Dy 1 * : protein ratio of 1:1. 
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(i) Assignment of pseudo-contact shifted peak of VI 7 : 

As shown in Fig. A 2 there are three new peaks (each marked by an asterix) 
compared to A], which arise due to PCS. It is impossible to assign these new peaks to 
any of the residues in £7zCaBP. For example, all these peaks serve as possible 
candidates for pseudo-contact shifted peaks of VI 7. This ambiguity is overcome by 
selective unlabeling of lie and Leu as shown in A 3 and A 4 . The spectral region is now 
simplified as selectively unlabeled amino acids don't show up their original and 
pseudo-contact shifted cross peaks. Thus, one is left with only cross peaks of VI 7, 
which can be unambiguously assigned to a lone pseudo-contact shifted cross peak 
seen in A 4 . Further, such spectral simplification helps in measuring the corresponding 
RDC (if any)as shown in A5. 

(ii) Assignment of pseudo-contact shifted peak of SI 33 (Fig. 2B) : 

As evident in Fig. 2B, simultaneous presence of cross peaks belonging to residues 
K28, K102, L38 and S133, result in ambiguity for the assignments of newly appeared 
pseudo-contact shifted peaks (each marked by an asterix in B 2 ). On selective 
unlabeling of Lys and Leu residues, cross peaks due to these residues are not seen in 
B 3 and B 4 . Thus, the new pseudo-contact shifted peak in B 4 can be unambiguously 
assigned to SI 33. The corresponding RDC for this residue is measured as shown in 
B 5 . 


(Hi) Assignment of pseudo-contact shifted peak of AU8 : 

In this case, two new peaks appear as a result of pseudo-contact shifted (each 
marked with an asterisk in C 2 ). However, from Ci and C 2 , it is not clear as to which of 
the new peaks correspond to the parent residue. However, removal of cross peaks 
belonging to LI 34 by selective unlabeling results in simplification of the spectra. 
From the earlier two cases shown in Fig. 2A and 2B, it is clear that new pseudo- 
contact peaks seem to appear with both 'H and 15 N chemical shifts up field shifted as 
compared to the parent cross peak. Thus, the new pseudo-contact shifted cross peak in 
C 4 can be unambiguously assigned to A 1 1 8 and the corresponding RDC is measured 
as shown in C5. 


Conclusion 

Although NMR spectroscopy is relatively a young technique compared to X-ray 
crystallography for structure determination of biomolecules, during the last decade it 
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has emerged as a powerful technique for structural studies of large proteins and 
protein complexes in solution. With significant recent improvements in techniques for 
over expressing isotope labeled proteins, which have taken place contemporaneously 
with the sophistication in spectrometer hardware and developments in new 
methodologies for improving spectral resolution and sensitivity, it is now possible to 
determine structures of bio-macromolecules more rapidly than ever before. In this 
direction amino acid selective "unlabeling" provides a useful methodology for 
structural studies of large proteins in solution. Significant simplification of the spectra 
is achieved using this approach, resulting in unambiguous resonance assignments. 
Further, this technique uniquely provides information in all the stages of structure 
determination of proteins by NMR. 
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Abstract 

The binding affinity of various ligands that bind to the histamine Hi receptor (PAT analogues) 
has been investigated using topological descriptors viz., negentropy (iV), molecular redundancy 
index (MRI), first-order molecular connectivity index (f), Wiener index ( W ), and Szeged 
index (Sz). The regression analysis gave a penta-parametric correlation involving N, MRI, Sz, 

IP! and IP 4 as the most significant QSAR model for modelling binding affinity of l-phenyl-3- 
amino-1, 2, 3, 4-tetrahydronaphthalenes as ligands. 

(Keywords : binding affinity / Hj receptor / QSAR / Wiener index / Szeged index / regression 
analysis) 


Introduction 

Quantitative structure-activity relationship (QSAR) is a major factor in 
contemporary drug designing. Thus, it is quite clear why a large number of users of 
QSAR are located in industrial research units. The coupling of QSAR methodology to 
the experience and intuition of professional drug designer should result in a much 
more organized search for the novel drugs of human, animal, and plant therapy 1-4 . 

The molecular descriptors like geometric, electronic, information theoretic, 
topological, and their combinations have widely been employed in QSAR studies 5-10 . 


# Paper presented at 70 th Annual Session of the National Academy of Sciences, India and 
dedicated to Vijay on his 43 rd birthday. 
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Booth et. al. n have shown that phenolamino tetralines (PATs) stimulate pyrosine 
hydroxylate activity and dopamine (DA) synthesis in rodent brain in-vivo. They have 
also shown 12 that this effect is medicated through interaction with a [N] binding site 
which was designated originally as the P AT-c> 3 -receptor, because although it has low 
affinity for G\/<s 2 ligands, if incorporated a c-like pharmacophore, it had a a-like 
rodent brain distribution, and its functional effects on dopamine synthesis were 
blocked by the BMY -14802. These data suggested a dependence between chirality 
receptor interaction, and functional effect. 

Recently, Bucholtz et at 3 have discussed that PATs also bind with high affinity to 
the [ 3 H]-mepyramine-labeled H] receptor in rat and guinea pig brain and also reported 
biological evaluation of additional PAT' antagonists (Table 1, Fig. 1) to identify 
differences in binding at these sites viz. [ 3 H] -(-)- trans-H 2 -PAT sites versus [ 3 H] 
mepyramine-labeled Hi a-receptors. However, til! date no quantitative structure- 
activity relationship (QSAR) was developed using topological indices. 

In view of the above, the present paper deals with topological modelling of 
binding affinity of l-phenyl-3-amino-I,2,3,4-tetrahydronaphthalenes as ligands for 
histamine H r receptors (Table 1, Fig. 1). 



Fig. 1 - 1 -phenyl-3 -amino-1, 2, 3, 4-tetrahydronaphthalenes used in the present study 
The binding affinity expressed as K 0 5 (nm) are taken from the literature 13. 
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Table 1- Structural details and activity of compounds used in the present study. 


Compd. 

No. 

Config. 

R. 

r 2 

r 3 

r 5 

R* 

r 7 

R^ 

Kos 

1 . 

(±)-trans 

H 

H 

N(CH 3 ) 2 

H 

Cl 

OH 

CH 2 

0.63 

2 . 

(±)-trans 

H 

H 

n(CH 3 ) 2 

H 

H 

H 

ch 2 

35.00 

3. 

(±)-trans 

H 

H 

NCH 3 (C 3 H 5 ) 

H 

H 

H 

ch 2 

3.40 

4. 

(±)-trans 

H 

H 

N(C 3 H 5 ) 2 

H 

H 

H 

ch 2 

10.20 

5. 

(±)-cis 

H 

N(CH 3 ) 2 

H 

H 

H 

H 

ch 2 

1200.00 

6 . 

(±)-trans 

H 

N(CH 3 ) 2 

H 

H 

H 

H 

ch 2 

940.00 

7. 

(±)-cis 

H 

H 

N(CH 3 ) 2 

H 

H 

H 

c 2 h 4 

109.00 

8 . 

(±}-trans 

H 

H 

N(CH 3 ) 2 

H 

H 

H 

c 2 h 4 

20.90 

9. 

(±)-trans 

H 

H 

N(ch 3 ) 2 

H 

OH 

OH 

ch 2 

60.00 

10 . 

(±)-cis 

H 

H 

N(CH 3 ) 2 

H 

OH 

OH 

ch 2 

9.60 

11 . 

(±)~trans 

ch 3 

H 

N(CH 3 ) 2 

H 

H 

H 

ch 2 

18.30 

12 . 

(±)-cis 

ch 3 

H 

N(CH 3 ) 2 

H 

H 

H 

ch 2 

1.90 

13. 

(±)-trans 

H 

H 

nch 3 

(CH 2 ) 2 C 6 H s 

H 

H 

H 

ch 2 

140.00 

14. 

(±)-trans 

H 

H 

nch 3 

(CH 2 ) 3 C 6 H 5 

H 

H 

H 

ch 2 

333.00 

15. 

(±)-trans 

H 

H 

nch 3 

(CH 2 ) 4 C 6 H s 

H 

H 

H 

ch 2 

212.00 

16. 

(±)-trans 

H 

H 

N(CH 3 ) 2 

o-Cl 

H 

H 

ch 2 

53.00 

17. 

(±)-trans 

H 

H 

N(CH 3 ) 2 

0 -CH 3 

H 

H 

ch 2 

10.20 

18. 

(±)~trans 

H 

H 

N(CH 3 ) 2 

p-Cl 

H 

H 

ch 2 

9.00 

19. 

(±)~trans 

H 

H 

N(CH 3 ) 2 

P-CH 3 

H 

H 

ch 2 

2.47 

20 . 

(±)-trans 

H 

H 

N(CH 3 ) 2 

p-F 

H 

H 

ch 2 

1.53 

21 . 

(±}-cis 

H 

H 

N(CH 3 ) 2 

H 

Cl 

OH 

ch 2 

0.54 


Contd.. 
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Table 1 Contd.. 


22. 

(±Hrans 

H 

H 

NH(CH 3 ) 

H 

Cl 

OH 

ch 2 

8.10 

23. 

(±)-trans 

H 

H 

NH(C 3 H s ) 

H 

H 

H 

ch 2 

45.00 

24. 

(±)-tr ans 

H 

H 

N(CH 2 ) 4 C 6 H s 

H 

H 

H 

ch 2 

1500.00 

25. 

(±)-trans 

H 

H 

nh 2 

H 

H 

H 

ch 2 

1270.00 


Methodology of Modelling Activity 

Pharmacology : 

In the present study we have adopted A!o .5 (nm) from the literature 13 . 

Non-linear regression analysis of inhibition data was used to determine IC50 
values using prism 2.0 (Graph Pat San Francisco, CA). IC 50 values were converted to 
K o s based upon the Cheng Prussof equation 14 , where K0.5 =IC 50 /(1 + L* / KD). 

Experiments were repeated at least three times to determine mean K 0 . 5 ± SEM. 

Molecular Descriptors : 


Molecular descriptors are calculated using the information given in the following 
Table: 


Molecular descriptors 

Abbreviation Formula (e) 

Basis 

Ref. 

Wiener index 

W= l/2±±d(v„v,p) 

M j=l 

Vertex distance in molecular 
graph 

8 

Szeged index 

sz= Y, n M G ) n M G ) 

e 

Vertex distance (near to edge 
ends) in molecular graph 

9,10 

Negentropy index 

N = (i).(n) where i = 
-K I Pi log?, 

j 

Information theory 

15 

Molecular 
redundancy index 

MRI- 

NlogN 

Information theory 

15 

First order valence 
connectively index 

ij 

Valence delta values 

1,2,16 
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The aforementioned topological indices are calculated from adjacency and 
distance matrices. The adjacency matrix, M, of a chemical structure is defined by the 
elements [My] where My is 1 if atoms / and j are bonded and zero otherwise. The 
distance matrix, D, of a chemical structure is defined by the elements [D tJ ] where Dy is 
the length of the shortest path from atom i to j ; zero is used if atoms i and j are not 
part of the same connected component. The adjacency matrix of CH 3 CKNO is 
displayed on the left and its distance matrix is displayed on the right (below). 



Structural graph of CH 3 CH=0 


Hydrogen suppressed molecular graph of CH 3 CH=0 


Ci 

H 2 

h 3 

H 4 


0 11110 0 
1 0 0 0 0 0 0 
1 0 0 0 0 0 0 
1 0 0 0 0 0 0 


0 1112 2 
1 0 2 2 3 3 

1 2 0 2 3 3 

1 2 2 2 3 3 - 


C 5 0 0 0 0 1 


0 


1 2 2 0 1 1 


H 6 


0 0 0 0 1 0 0 


2 3 3 1 0 2 


0 7 


1 0 0 0 1 0 0 


2 3 3 1 2 0 


The adjacency matrix and distance matrices of chemical structures are promising 
techniques in molecular modelling. In graph theory the molecular structure is 
transformed into molecular graph by depleting all the carbon and hetero-atoms 
hydrogen bonds. Hence, hydrogen suppressed molecular structure is defined as 
molecular graph. 


A number of successful QSAR studies were made based on the Wiener index, 
however, very little is known on the use of Szeged index in developing quantitative 
structure-activity relationships. 



288 


VIJAY K. AGRAWAL et. al. 


Table 2- Calculated molecular descriptors and indicator parameters. 


Compd. 

No. 

N 

y 

MRI 

Sz 

W 

ip] 

Ip 2 

Ip 3 

Ip 4 

1. 

52.9074 

0.224 

8.2016 

1960 

1072 

1 

0 

1 

1 

2. 

50.7645 

0.3176 

7.7072 

1524 

826 

I 

0 

0 

1 

3. 

58.5225 

0.2133 

8.0286 

1692 

936 

0 

0 

0 

1 

4. 

64.4987 

0.2212 

8.6201 

1911 

1139 

0 

0 

0 

1 

5. 

47.2402 

0.2856 

7.1667 

1466 

730 

0 

0 

0 

0 

6. 

47.2402 

0.2856 

7.1667 

1328 

718 

0 

0 

0 

1 

7. 

50.1908 

0.3059 

7.6791 

1434 

798 

1 

0 

0 

0 

8. 

52.5976 

0.2726 

7.6892 

1536 

834 

1 

0 

0 

1 

9. 

53.1609 

0.2431 

6.7013 

2183 

1207 

1 

0 

1 

1 

10. 

50.7486 

0.2771 

7.4413 

2159 

1207 

1 

0 

1 

. 0 

11. 

51.9948 

0.281 

7.5446 

1490 

830 

1 

o 

0 

1 

12. 

51.9992 

0.2809 

7.5874 

1360 

792 

1 

0 

0 

0 

13. 

70.308 

0.2484 

9.81 

2776 

1852 

0 

0 

0 

1 

14. 

76.2432 

0.2382 

10.31 

3454 

2125 

0 

0 

0 

1 

15. 

80.442 

0.246 

10.81 

3822 

2426 

0 

0 

0 

1 

16. 

49.4132 

0.2527 

7.6952 

1531 

829 

1 

1 

0 

1 

17. 

54.1728 

0.2502 

7.5933 

1531 

829 

• 1 

1 

0 

1 

18. 

49.4132 

02527 

7.6893 

1610 

867 

1 

1 

0 

1 

19. 

54.1728 

0.2502 

7.5875 

1610 

867 

1 

1 

0 

1 

20. 

49.4132 

0.2527 

6.9756 

1610 

867 

1 

1 

0 

1 

21. 

50.5008 

0.2593 

7.8196 

1942 

1072 

1 

0 

1 

0 

22. 

50.6181 

0.1843 

7.8377 

1960 

1072 

0 

0 

1 

1 

23. 

53.1342 

0.2206 

7.5122 

1732 

976 

0 

0 

0 

1 

24. 

69.2493 

0.3081 

10.4306 

3822 

2426 

0 

0 

0 

1 

25. 

39.823 

0.2631 

6.3574 

1274 

733 

0 

0 

0 

1 


TOPOLOGICAL MODELLING OF BINDING AFFINITY OF SOME LIGANDS 


289 


Several computer programs are available for computing molecular descriptors 
including topological indices. Some of the popular packages in calculating molecular 
descriptors are: AHPAL, CODESSA, MOLCONN-X, HYPERCHM, DRAGON, 
POLLY, MOPAC, ALCHEMY III, TSAR, HASL, CAChe, Chem-3D, TOPMOST 
CLOGP. In addition, several softwares are available to carry out regression analysis: 
MSTAT, SPSS-X, BMD02R, BMLR, SPSS/PC, MINITAB, SAS etc. 


Indicator parameters (Ip u Ip 2, IP 3, Ip 4) •' 

Indicator variables 17,26 (parameters), are some times used in linear multiple 
regression analysis to account for certain features which cannot be described by 
continuous variables. In QSAR equations they normally describe a certain structural 
element, be it a substituent or another molecular fragment. Thus, Free Wilson analysis 
may be interpreted as a regression analysis approach using only indicator variables. 

The disadvantage of indicator variables are that they are dummy variables and that 
too many of them increases the possibility of chance correlation as well as collinearity 
defect 18,21 - 26 . 

In the present case the indicator parameter Ipi was taken as unity when -N(CH 3 ) 2 
group is present at R 3 . Similarly, Ip 2 is used as unity when R 5 is other than H. When 
R^ and R 7 are other than H the indicator parameter taken isIp - 3 and its value is taken 
as unity. Ip 4 is taken as one when the compound is trans. 


Regression Analysis : 

We have used the maximum R 2 improvement method to identify prediction 
models 17-26 . This method finds the "best" one variable model, the "best" two variable 
model and so forth for the prediction of property/activity. Several models 
(combinations of variables) were examined to identify combinations of variables with 
good prediction capabilities. In all regression models developed we have examined a 
number of statistical parameters associated with residues, i.e. the Wilks-Shapiro test 
for normality and Cooks D-statistics for outliers, to obtain the most reliable results 26 . 

Multiple regression analysis for correlating binding affinity (X 0.5 Value) of the 
present set of compounds with the aforementioned molecular descriptors were carried 
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out using Regress-1 software as supplied by Professor I. Lukovits, Hungarian 
Academy of Sciences, Budapest, Hungary. Several multiple regressions were 
attempted using correlation matrix from this program and the best results were 
considered and discussed in developing QSAR. 

Computations : 

All the computations were carried out in Power Macintosh 9600/233. 

In the present study, we have used Wiener-1 and Regress-1 programs for the 
calculation of W and Sz indices and for making regression analysis respectively. Other 
calculations are made for programmable calculators. These softwares were made 
available to us by Prof. I. Lukovits, Hungarian Academy of Sciences, Budapest, 
Hungary. 


Results and Discussion 

A perusal of Table 1 shows that no correlation exists between configuration and 
the inhibition constant (binding affinity constant), Kq.s (nm). That is, we can not use 
this parameter for distinguishing between cis-and trans-configurations. Surprisingly 
enough, we observed that the topological indices used (N, MRI, 'x K , Sz, and W) have 
different values for these conformations (Table 1). It is observed that the values of 
these topological indices are higher for trans-conformers compared to the cis- 
conformers. Hence, these indices are appropriate for modelling binding affinity of 1- 
phenyl-3-amino-l,2,3,4-tetrahydronaphthalenes as ligands for histamine Hi -receptors. 

In order to obtain statistically significant models, we have employed the multiple 
regression analysis technique and obtained the correlation matrix as shown in Table 3. 
The descriptors in the selected sub set are those which exhibit correlation coefficient 
below 0.5. 

A perusal of Table 3 shows the existence of collinearity between: (i) N and '% v (ii) 
N and Sz; (iii) N and W; (iv) V and Sz; (v) V and W; (vi) Sz and W. The QSAR 
models based on such combinations of topological indices will suffer from the defect 
of collinearity. None of the indicator parameters (Ip 1; Ip 2 , Ip 3 and Ip 4 ) correlate 
mutually or with any of the topological indices (N, MRI , Sz and W) used. 
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Table 3- Correlation matrix for the intercorrelation of molecular description and correlation with the 
activity. 



K 0 5 

N 

y 

MRI 

Sz 

W 

ip. 

ip 2 

ipj 

Ip4 

K os 

1.0000 










N 

-0.0202 

1.0000 









MRI 

0.4013 

-0.1811 

1.0000 








V 

0.0976 

0.9447 

-0.0737 

1.0000 







Sz 

0.2097 

0.8836 

-0.0862 

0.8974 

1.0000 






w 

0.2232 

0.8990 

-0.0654 

0.9149 

0.9948 

1.0000 





Ip. 

-0.5478 

-0.4335 

0.3041 

-0.4358 

-0.4228 

-0.4486 

1.0000 




Ip 2 

-0.2523 

-0.2035 

-0.0917 

-0.2185 

-0.2556 

-0.2679 

0.4432 

1.0000 



Ipj 

-0.2517 

-0.1892 

-0.3203 

-01775 

0.0636 

0.0048 

0.2418 

-0.2500 

1.0000 


Ip 4 

-0.0275 

0.2662 

-0.3941 

0.2049 

0.1908 

0.2002 

-0.2418 

0.2500 

-0.2500 

1.0000 


The values of correlation coefficient in Table 3 also show that none of the 
molecular descriptors correlate with the inhibition constant (affinity constant) K 0 .s 
(nm). However, marked relationship occurred between MRI, Ipi with the K 0.5 nm. 
This indicates that statistically significant multi-parametric model will contain two 
molecular descriptors. 

Applying multiple linear regression technique we have obtained several 
regression models. Only those models wherein the standard deviation of the 
coefficient is smaller than the values of coefficient of the descriptors are considered as 
statistically significant models. Such models are only considered and discussed 
critically using the data presented in Table 4. 


The data presented in Table 4 shows that : 
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Table 4- Quality of correlation and other statistical paparameters of significant correlations. 


Model 

No. 

Regression equation 

SE 

R 

F-ratio 

Q 

1. 

*o. 5 = (-14.99) AH614.23) Ip,+1410.24 

372.85 

0.6180 

6.796 

0.0017 

2. 

*o. s = (8998.06) MRI-(661.50) Ip, -1705.69 

278.37 

0.8096 

20.924 

0.0029 

3. 

*05 = (-13.24) ^+(8779.16) MRI-(766.21) Ip,-860.74 

257.31 

0.8479 

17.909 

0.0033 

4. 

*o.5 = (-87.52)V + (9203.65)MRI-(75 1 .59) Ip, -1008.1 1 

268.20 

0.8335 

15.928 

0.0031 

5. 

*o. 5 =(9633.23)MRI-(743.83)Ip, + (196.81)Ip2-1862.43 

275.08 

0.8239 

14.795 

0.0030 

6. 

*o.s=U 0027. 19)MRl-(7 1 6.73)Ip,+( 1 83. 42)Ip 3 -l 976.32 

276.12 

0.8224 

14.631 

0.0030 

7. 

*o. s = (-42.33) N - (7822.09) MRI + (0.44) Sz - 
(718.47) Ip, + 114.35 

208.46 

0.9079 

23.463 

0.0044 

8. 

* 0.5 — (—44.60) N + (7406.58) MRI + (0.67) VV - 
(693.09) Ip, + 428.79 

211.78 

0.9048 

22.579 

0.0043 

9, 

ATo 5 = (-13.00) N + (9386.81) MRI - (842.62) Ip, + 
(187.09) IP 2 - 1024.61 

253.56 

0.8603 

14.238 

0.0034 

10. 

^o. 5 = (—42.33) N- (7822.09) MRI + (0.44) Sz- 
(718.47) Ip, -35.68 

193.99 

0.9249 

22.494 

0.0048 

11. 

*o.5 = (-57.41) N+ (5788.88) MRI + (0.64) Sz- 
(624.63) Ip, - (278.45) Ip 3 + 1077.60 

198.32 

0.9214 

23.259 

0.0047 

12. 

*o.5 = (-43.85) N + (8399.77) MRI + (0.44) Sz - 
(71 1.16) Ip, + (133.77) Ip 4 - 76.54 

206.71 

0.9143 

19.359 

0.0044 


(1) No mono-parametric models are possible for modelling the activity (K 0 5 
nm); 

(2) Only two bi-parametric models are obtained and the one using MRI and Ip, 
gave better value for R (= 0.8096); 
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(3) All the four tri-parametric models are statistically significant and the 
triparametric correlations containing MRI, N and Ip, gave the highest value 
ofR(= 0.8479). 

Comparing this result with the bi-parametric correlations discussed above we 
observed that the introduction of N to MRI and Ip, resulted into improvement of 
quality of model such that R value changed from 0.8096 (bi-parametric) to 0.8479 (tri- 
parametric). 

(4) Three tetra-parametric models are possible which are statistically significant 
and each gave improved R-values compared to the best triparametric model 
discussed above. The tetra-parametric models; (i)involving MRI, N, Sz and 
Ip, and (ii) involving MRI, N, V/. and Ip, have similar quality as their R- 
values are identical (0.9079 and 0.9048). These identical . qualities may be 
attributed to strong inter-correlation between W and Sz. Because of this 
information content in W and Sz is same. 

(5) Finally, step-wise regression resulted into three penta-parametric 
correlations. All these three models are statistically excellent. However, the 
model containing MRI, N, Sz, Ip, and Ip 2 has the highest value of i?.Second 
best penta-parametric model is the one containing MRI, N, Sz, Ip, 'and Ip 3 , 
(i? = 0.9214). It is worthy to record that none of these models contained W as 
one of correlating parameters. This indicates that Sz is a better topological 
index than W in proposing multi-parametric models. 

It is interesting to record that all the proposed models contain invariably MRI and 
Ip, as one of the correlating parameters. Also that the coefficient of MRI is always 
positive and that of Ip, is always negative. Recall that Ip, is used for the presence or 
absence of -N(CH) 2 moiety at R 3 . The results, therefore, indicate negative effect of 
this group on the activity (Kq.s, nm) of the compounds used in the present study. 


In order to confirm our aforementioned findings we have calculated quality factor 
g 27 . This quality factor, Q, is defined as the ratio of the correlation coefficient R to the 
standard error of estimation, SE (Q = ft/SE).That is, now the quality of model is 
judged considering R and SE simultaneously. Hence, Q is a better parameter than R or 
SE to decide the quality of models. Such Q values are given in Table 4. 
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The 0-values presented in Table 4 indicate that the quality of models (7), ( 8 ) (10) 
(11) and (12) are similar. Hence, in the present case 0-value failed to give which 
model is most suitable for modelling the inhibitor binding affinity constant. Hence, 
one has to use cross-validation method for deciding predictive ability of the models. 

Another alternative for obtaining the most suitable model is to estimate K 0 . 5 and 
compare them with its observed values and finally obtain predictive correlation 
coefficient (R 2 ) .by the regression analysis of correlation between observed and 
estimated K Q 5 from different proposed models. 



Fig. 2- Correlation between observed and estimated activity using various proposed models. 


We have, therefore, calculated K 0.5 values using models ( 7 ), ( 8 ), (10), (11) and 
(12) and obtained R 1 values, highest R 2 ( 0.986) is observed for the model (12) 
indicating it to be the most appropriate model for modelling binding affinity (K o s , nm) 
of the series of compounds used in the. present study. 
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Abstract 

A series of the complexes of the type [Co(DH) 2 (N 3 )L] where DH = dimethyl glyoxime and L = 
pyridine or its derivatives have been synthesized and characterized on the basis of electronic, 

IR, NMR ('H and 13 C) and FAB MS as well as thermal analysis. IR and NMR spectra indicate 
the trans configuration for the isolated complexes. The FAB MS shows molecular ion peak 
which corresponds to mono nuclear complex of the type [Co(DH) 2 N 3 L]. 

(Keywords : azido cobaloximes/ dimethylglyoxime/ pyridine/ spectral studies/ trans 
configuration) 

Introduction 

The bis(dimethylglyoximato)cobaIt(III) complexes, trivially known as 
cobaloximes in view of their resemblance to cobalamins. There has been growing 
interest in organo cobalt complexes which might serve as ‘model compounds’ for 
coenzyme Bn These have been extensively studied in the last three decades 1 ' 9 in view 
of their rich chemistry, both with respect to synthesis as well as catalysis 10 ' 14 . 
Synthesis of organo cobaloximes with modified structural features continues to be a 
fascinating area in the chemistry of model compounds of coenzyme B ]2 . Recently, 
Gupta and Roy 15 reviewed the synthetic aspects of organo cobaloxime chemistry. 
Further, the X-ray crystal structure of alkyl cobaloximes suggests 16 that the Co atom 
in cobaloximes has a slightly greater positive charge, causing relatively stronger 
attachment of a base along axial directions.. In addition to this the azido ligand present 
in cobaloximes shows many interesting chemical features 17,18 . All these 
considerations prompted us to synthesize a series of azido cobaloximes containing 
pyridine and substituted pyridines. We report here the synthesis, IR, electronic, 
NMR( ! H and 13 C), mass spectral data as well as thermal analysis of cobaloximes 



298 


N. RAVI KUMAR REDDY et. al. 


containing pyridines as axial ligands and the linear tri-atomic azide ion (Nf) in the 
sixth coordination position. 


Materials and Method 

Cobalt(II)acetate , KN 3 , dimethylglyoxime were obtained from BDH . Pyridine 
and substituted pyridines (4-aminopyridine, 4-picoline, 4-ethylpyridine, 4-bromo 
pyridine hydrochloride, 4-cyanopyridine, 4-vinylpyridine and isonicotinamide) were 
used as supplied by SD Fine chemicals. 

Synthesis of [Co(DH) 2 (N 3 )L] complexes : Complexes of the type 

[Co(DH) 2 (N 3 )L] (where L = pyridine and substituted pyridines) were prepared in 60- 
70% yield by bubbling air for a period varying from 2 to 6 h through an ethanolic 
solution of cobalt(II)acetate, KN 3 , dimethylglyoxime and the ligand taken in a 1: 2: 2: 
2 proportion. This was allowed to stand for to settle the micro crystalline solid at the 
bottom and collected on a frit, washed with small aliquots of EtOH, Et 2 0 and finally 
dried in a vacuum dessicator. The reaction presumably proceeds in the following 
manner. 


2Co(CH 3 COO) 2 .4H 2 0 + 4MN 3 + 4DH 2 + 4L + 'A 0 2 ► 

2[Co(DH) 2 (N 3 )L] + 4CH 3 COOM + 2L.HN 3 + 9H 2 0 

Synthesis of [Co(DH)2(N 3 )(4-NH 2 Py] : 0.5g (2mmol) Co(OAc) 2 .4H 2 0 was 
dissolved in 25ml ethanol. 0.33g (4.1mmol) KN 3 was added to the above solution with 
constant stirring. Then 0.5g (4.3mmol) dimethylglyoxime was added to the mixture 
and heating and stirring continued for 30min till a deep brown solution was obtained. 
To the filtrate 0.395g (4.2mmol) of amino pyridine solution (dissolved in 5ml ethanol) 
was added and stirred for 5 more min and the solution allowed to cool at room 
temperature. A continuous stream of air was then bubbled through the solution for 2'A 
h when the reddish-brown micro crystalline solid was deposited at the bottom. After 
allowing the precipitate to settle, it was collected on a frit, washed with ethanol 
several times followed by dry ether. The product (73%) was finally dried in vacuo. 

Cobalt was determined by the method of Laitenen and Burdett 19 . Micro analysis 
of complexes were performed on Perkin-Elmer 240C elemental analyzer. IR spectra 
were obtained on Perkin Elmer FTIR 1605 instrument using KBr pellets. 'H and l3 C 
NMR spectra were recorded on Varian Gemini 200 NMR spectrometer and referenced 
to TMS in CDCl 3 /DMSO-d6. Mass spectra were recorded on Autospec mass 
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spectrometer and UV -Visible spectra of the complexes were recorded on Hitachi U- 
3410 spectrophotometer. Thermal analysis was obtained on Mettler Toledo Star 
system. 


Results and Discussion 

All these complexes have very low molar conductivity values which favour the 
neutral, non-ionic formulations and non-electrolytes in dilute solutions. The ligand- 
field spectra of the complexes in MeOH show a peak of w'eak to moderate intensity at 
around 18000-20000 cm' 1 20 . There are two other main bands centered at 40000— 
42500 cm -1 region attributable to the intra-ligand n-n* transition of the coordinated 
oxime group, where as the band at 28000-30000 cm' 1 may be due to the pyridine — >• 
cobalt transition. The strong band at around 47000 cm -1 due to the d7i(Co)-rc*(DH“) 
transition 21,22 . The general features in the spectrum of these complexes are very 
similar to those observed for the corrins and other alkyl cobaloximes 23 . 

The IR spectra of the free ligands and the complexes have been carefully analysed 
since this technique serves as a probe to distinguish between the cis- and trans- forms 
in organo cobalt(III) compounds 24,25 and these are all shown to have indeed the trans 
configuration similar to the organo cobaloximes 26 , which is further substantiated on 
the basis of NMR evidences. Most of the important bands of all these complexes 
appear in almost similar positions with identical intensity depicting the similarity in 
structure. The free uncoordinated DH 2 has no band at 1240 cm"', a strong band 
appears in the Ni(DH) 2 complex at 1240 cm' 1 due to u(N-0) of the ionised N-OH 
group of dimethylglyoxime 27 . In the present case, for all azido cobaloxime complexes, 
presence of two bands in the region 1230-1240 cm' 1 and 1090-1095 cm -1 account for 
two non-identical linkages i.e. C=N-0-H and C=N-0 H in the complexes 4 . This 
observation is consistent with the X-ray diffractional analysis of the Cu(DH) 2 
complexes in which two different Cu-N distances were observed 28,29 . The stretching 
vibration, u(C=N), of the non-coordinated DH 2 appearing at 1620-1 630cm" 1 is 
consequently shifted in this azido cobaloximes to 1560-1570 cm -1 . A very weak 
intensity band observed in the region 1700-1750 cm' 1 is attributed to the deformation 
of the strong intramolecular hydrogen bridge as it disappears on deuteration 30,31 . The 
characteristic u(Co-N) (N of DH~) appears at 510-515 cm' 1 and u(Co-N)(N of 
pyridine) at 420 cm' 1 . The IR bands of the N 3 ' ion , which is a linear symmetric tri 
atomic system with sixteen valency electrons, are of interest. The bending and 
strectching modes of the azide group 32,33 appearing at 640 and 2040 cm -1 are shifted to 
660-680 cm -1 and 2015-2025 cm' 1 , respectively. N 3 " stretching frequency confirms 
that the complex is mono nuclear and N 3 ' coordinated to Co(III) in end on mode 34 . 
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The characterstic absorption bands in dimethyl glyoximato ligands in these 
cobaloximes are also sensitive to the change of axial ligands. As the donating power 
of the base ligand increases the uO-H...O at 1770 cm -1 and oC=N at 1560 cm 1 shift 
to lower wave number region while o(N-0) at 1230 cm 1 and 1085 cm 1 shift to 
higher one. These shifts can be interpreted as follows. The coordination of more 
electron donating base to Co(III) causes the increase in electron density in Co(III) 
which facilitates the back donation of Co(III) to the nitrogen atom of dimethyl 
glyoximato ligand, resulting in the increase in electron density in C=N and N-O 
bonds. The increase in electron density in N-0 bonds causes the stronger hydrogen 
bridges 0-H...0 and higher frequency shift of N-0 stretching vibration. The 
facilitated back donation from cobalt(III) to nitrogen atoms of dimethyl glyoxime 
means the increased metal donor n bond in the equatorial moiety of cobaloximes, 
which causes the stronger interaction of Co with equatorial ‘N’ atoms resulting in 
higher frequency shifts of o(Co-N) vibration and causes more conjugation in the five 
membered chelate ring including Co(III) which effects the u(C=N) vibrations to shift 
to lower frequency. 

The NMR (’H and 13 C) spectra have also been used in elucidating the stereo 
chemistry and in assessing the interaction of ligands coordinated to the same metal 
ion 35 . The 'H NMR spectra of the azido cobaloxime complexes clearly show that on 
coordination, the signal due to CH 3 group of DH 2 at 1.98 ppm undergoes a clear down 
field shift by around 0.40 ppm to 2.36-2.42 ppm. In all these complexes, this signal 
appears as a singlet indicating the equivalence of four methyl groups and presence of 
trans structure 9 . The NMR spectra of simple free ligand of pyridine and other 
derivatives have been assigned. The aromatic protons of pyridine compounds absorb 
in the aromatic region. In free pyridine the C 4 -H proton to the hetero atom observed at 
7.75 ppm as triplet, the C 2 -H(C 6 -H) protons and C 3 -H(Cr-H) protons observed at 
7.60 ppm and 7.25 ppm respectively as multiplets. On coordination, pyridine to 
Co(III) in metal complex, these protons signals observed at 7.82 ppm, 8.18 ppm and 
7.39 ppm for C 4 -H, C 2 -H(C 6 -H) and C3-H(C5~H) respectively. The more down field 
shift of C 2 -H(C 6 -H) protons clearly shows the coordination of N atom of pyridine to 
the Co atom. 4-Bromo pyridine free ligand proton chemical shifts observed at 7.80 
ppm and 7.30 ppm for Cr-H(Cs-H) and C 3 -H(Cr-H) protons, respectively. On 
complexation, these protons undergo down field shift and observed at 8.16 ppm and 
7 .42 ppm, confirms the formation of the complex through N of pyridine ring. In all the 
other complexes the chemical shifts follow the same trend and the data are presented 
in Table 3. 



Table 1- Analytical data of [Co(DH) 2 (N 3 )L] complexes 
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Table 4- l3 C ['H]NMR spectral data* of [Co(DH) 2 (N 3 )L] complexes 
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13 C NMR spectroscopy has been recognised as one of the most promising tools 
for the study of vitamin-B^ and its model compounds 36 ' 39 . A single peak was 
observed as expected for azido cobaloximes at 12.70—12.90 ppm and the imine carbon 
chemical shift at ~150.0 ppm. This confirms the trans configuration for these 
complexes. The 13 C signals of the free pyridine ligands undergo downfield shift on 
complexation. The C 2 (C 6 ),C 3 (C 5 ) and C 4 signals of pyridine ring resonating at 148.9 
ppm, 124.5 ppm and 135.7 ppm respectively. In the [Co(DH) 2 (N 3 )(Pyridine)] 
complex, these values are shifted downfield and observed at 152.62 ppm, 128.35 ppm 
and 140.62 ppm respectively. In free 4-ethyl pyridine, the CH 2 and CH 3 signals are 
observed at 27.10 ppm and 12.70 ppm. These signals in the complex appeared at 
27.70 ppm and 13.35 ppm with slight down field. All these complexes presented in 
this work follow the same trend of down field shift. 13 C NMR chemical shifts for all 
the complexes studied are reported in Table 4. 

The mass spectral studies of free DH 2 , Ni(DH )2 complexes reveal the formation of 
molecular ion peaks for both DH 2 and Ni(DH) 2 40,41 . The FAB MS studies of 
[Co(DH) 2 (N 3 )(4-BrPy)] reveal the formation of molecular ion peak at 490. This 
corresponds to M +1 of [Co(DH) 2 (N 3 )(4-BrPy)]. Peak at 289 corresponds to the 
Co(DH) 2 moiety. FAB MS supports that complex is mono nuclear. The thermal 
decomposition studies indicate that thermal stability of all these cobaloximes is very 
high and with increase in temperature mass losses are seen corresponding to 
polymeric intermediates and the ultimate formation of Co 3 C >4 as the end product 40-42 . 
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Abstract 

Four new complexes of bivalent transition metals Co, Ni, Zn and Cd with schiff base, 2- 
hydroxy-5-methyl benzene-l-carbalidene-2-aminothiazoIe-3-carbalidene-2-aminobenzothiazole 
(HMBCATCABT) having the stoichiometry ML 2 , where M = Co (II), Ni (II), Zn (II) and Cd 
(II) and L= ligand have been synthesized. They have been characterized by their repeated 
melting point determination, T.L.C. for single spot, elemental analyses, conductance 
measurements, magnetic susceptibility, IR, NMR and UV/VIS spectral data. 

(Keywords : heterocyclic schiff base/IR/NMR/UV/VIS spectra). 

Introduction 

The chemistry of metal complexes with schiff base ligands has gained much 
importance due to their diversified biological applications 1 * * ' 5 . They are also found to 
have a number of pharmacological utilities 6 . This fact prompted us to synthesize some 
new schiff bases and their metal complexes to get better understanding in this field. 

In continuation of our earlier work 7 , we have prepared and characterized four 
complexes of different transition metals with the synthesized ligand (2-hydroxy-5- 
methyl benzene- 1 -carbalidene-2-aminothiazole-3-carbalidene-2-aminobenzothiazole). 

Materials and Method 

All chemicals used were of A.R. grade and the solvents were purified and distilled 
prior to their use. 

1. Synthesis of 2-hydroxy-5-methylbenzene-l ,3- dicarbaldehyde (HMBD) : 2- 

Hydroxy-5-methylbenzene-l,3- dicarbaldehyde was prepared according to the 

procedure given by Gagne et al. z [Yield-4.5 g (28%)]. 
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2. Synthesis of 2-hydroxy-5-methylbenzene-l-carbalidene-2-amino thiazole-3- 
carbalidene-2-aminobenzothazole (HMBCA TCABT) .-To 15 ml methanolic solution 
of 0.16 g (ImM) of dicarbaldehyde, a mixture of 0.10 g (ImM) of 2-aminothiazole 
and 0.15 g (ImM) of 2-aminobenzothiazole in 30 mimethano! was added. The 
resulting mixture was refluxed for 8-10 h. The excess of methanol was distilled off. 
On cooling, bright orange crystalline solid was isolated which was washed with 
alcohol, recrystallized from methanol and dried in vacuum desiccator over anhydrous 
CaCl 2 [Yield- 0.25 g (67%)]. 

3. Synthesis of metal complexes : To 20 ml methanolic solution of ligand 
(HMBCATCABT, 0.76 g, 2 mM), 15 ml aqueous solution of 1 mM of metal salt [0.25 
g,(CH 3 COO) 2 Co.4H 2 O/0.25g, (CH 3 COO) 2 Ni.4H 2 O/0.22 g, (CH 3 COO) 2 Zn. 2 H 2 O/0.18 
g, CdCl 2 ] was added and refluxed for 3 h. On cooling, coloured solid separated out. 
The resulting complex was filtered, washed with alcohol, ether and recrystallized from 
DMF. It was dried in vacuum desiccator over anhydrous CaCl 2 

Physical and analytical measurements : Carbon, Hydrogen and Nitrogen analyses 
were carried out on Carlo Erba Micro Analyser (Model- 1106) and Sulphur was 
estimated as BaS04 by standard procedure 9 . Metals were estimated by precipitating 
them as their pyridine complex 10 . Molecular weights of all the compounds were 
determined by Cryoscopic method 11 . IR spectra were recorded on a JASCO 
Spectrometer (Model IR-Report-100) in KBr medium. r H NMR spectra were carried 
out on NMR Brucker AC-300 MHz in CDC1 3 . The UV/VIS spectra in DMSO were 
recorded on Shimadzu Digital Double Beam Spectrophotometer (Model UV 150- 
150.02). The conductance measurements were carried out on Toshniwal Conductivity 
Bridge using 10 3 M DMF solution. The magnetic susceptibility was determined at 
room temperature by Gouy's method. 


Results and Discussion 

The analytical data of complexes reveal their 1 :2, (ML 2 ) stoichiometry and the 
conductivity measurements suggest their non-electrolytic nature. 

In IR spectra of dicarbaldehyde, sharp bands in 1715-1695 cm' 1 and at 2880 cm" 1 
were obtained due to the stretching vibrations of carobnyl 15 and -CH-(aldehyde) 13 
groups respectively. The bands due to phenolic group 14 and -CH-(aromatic) were 
obtained at 3525 cm -1 and 3000 cm" 1 respectively. 



Table 1-Physical and analytical data. 
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The infra-red spectra of schiff base (HMBCATCABT) show two bands in the 
region of 695-675 cm' 1 and 650-630 cm' 1 which may be due to symmetric and 
asymmetric vibrations of C-S stretching of thiazole or benzothiazole 15 while the C = N 
(cyclic) band has been observed in the region of 1485-1460 cm' 1 . The sharp intensity 
bands at 1625 cm' 1 and 3530 cm' 1 due to stretching vibrations of azomethine group 16 
and phenolic group were obtained. In the 1R spectra of metal complexes, the shifting 
of azomethine band indicates the coordination of ligand to metal through nitrogen. 
The band due to phenolic group was absent in the complexes indicating the 
deprotonation of phenolic group. However, no shift in the band position due to 
sulphur (thiazole and benzothiazole ring) was observed which infers that the sulphur 
of the ligand does not coordinate. 


1 H NMR spectra in CDC1 3 , of the synthesized dicarbaldehyde show 8 values at 2.2 
(3H,s Me), 7.8 (2H,m, Ar-H) 17 ’ 18 , 10.2 (2H,s, aldehydic-H) 19 and 11.4 (1 H,s, 
phenolic-H). In NMR spectra of ligand (HMBCATCABT), new 5 value at 8.5 (2 H,s, 
azomethine-H) 20 was observed showing the condensation of dicarbaldehyde with 
ligand. A sharp signal characteristic of phenolic proton in ’H NMR spectra of free 
ligand was found to be absent in their corresponding metal complexes. This 
observation confirms the deprotonation of phenolic group in the metal complexes and 
supports the formation of M-0 bond. 


In the electronic spectra of Co (II) complex, three bands were obtained at 8095 
cm' 1 , 15785 cm' 1 and 18945 cm' 1 corresponding to transitions 4 7i g — » 4 7V-> 
^ T 1 g - T ig ^ A 2 g and 4 r !g (P) respectively. The 10 Dq and P values were found at 
7690 cm' 1 and 696.3 cm' 1 respectively. The UV/VIS data suggests an octahedral 
geometry 21 around the metal ion. The magnetic moment value at 4.14 B.M. 22 of the 
complex also supports the octahedral geometry. The appearance of three bands in the 
electronic spectra of Ni (II) complex at 11130 cm" 1 , 17115 cm" 1 and 26265 cm' 1 
correspond to transitions 3 ^ 2 g-> 3 P 2 g, 3 ^ 2 g-> 3 Pi g and 3 7\ g (P) respectively. These 
transitions suggest an octahedral geometry 2 around the central metal ion which is also 
supported by its effective magnetic moment value al 3.18 B.M. The 10 Dq and p 
values were obtained at 11130 cm' 1 and 666 cm' 1 respectively. 
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Proposed structure of ligand (HMBCATCABT) 

ch 3 



ch 3 

[Where M 2+ = Co 2+ , Ni 2 % Zn 2+ and Cd 2, [ 


Proposed structure of metal complex 
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Abstract 

The propagation of surface waves in a cylindrical cavity embedded in a micropolar generalized 
thermoelastic medium is discussed in two cases (i) the cylindrical bore is empty (ii) the 
cylindrical bore is filled with homogeneous inviscid liquid. The dispersion equation appropriate 
to the surface wave propagation near the surface of the cylindrical bore is obtained, 
characterizing the dispersive nature of _ the wave. Some special cases have been deduced. The 
numerical results obtained have been illustrated graphically to understand the behaviour of 
phase velocity versus wave number of the wave. Significant thermal and microrotationa! effects 
are observed. 

(Keywords : surface wave propagation/micropolar/ generalized thermoelasticity) 

Introduction 

The problem of surface wave propagation near the bore hole in elastic host 
medium is of great practical importance. Since valuable organic and inorganic 
deposits beneath the earth's surface are difficult to detect by drilling randomly, the 
wave propagation technique is the simplest and most economic and does not require 
any drilling through the earth. Thermoelasticity theories which admit a finite speed for 
thermal signals (second sound) have aroused much interest in the last three decades. 
Recently relevant theoretical developments in the theory of thermoelasticity on the 
subject of .finite velocity of heat propagation are due to Green an Naghde 1 ' 3 , which 
provide sufficient basic modifications in the constitutive equations that permit 
treatment of a much wider class of heat flow problems. An important feature of 
this theory, which is not present in other thermoelasticity theories is that this 
theory does not accommodate dissipation of thermal energy. 

"Micropolar elasticity" termed by Eringen 4 is used to describe deformation of 
elastic media with oriented particles. A micropolar continuum is a collection of 
interconnected particles in the form of small rigid bodies undergoing both 
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translational and rotational motions. Typical examples of such materials are granular 
media and multimolecular bodies, whose microstructures act as an evident part in their 
macroscopic responses. The physical nature of these materials needs an asymmetric 
description of deformation, while theories for classical continua fail to accurately 
predict their physical and mechanical behaviour. For this reason, micropolar theories 
were developed by Eringen 4 ' 15 for elastic solids, fluids and further for non-local polar 
fields and are now universally accepted. 

The linear theory of micropolar thermoelasticity was developed by extending the 
theory of micropolar continua to include thermal effects by Eringen and Nowacki 8 . 
Different authors 9 " 14 discussed different type of problems in elasticity/ 
micropolar elasticity/micropolar generalized thermoelasticity. The present 
study is concerned with the problem of surface wave propagation in a 
micropolar generalized thermoelastic solid containing a cylindrical bore. 

Problem Formulation and Solution 

Case I-Propagation of Waves in cylindrical empty bore : Consider a cylindrical 
empty bore of radius 'a' having a circular cross-section in a homogeneous, isotropic 
micropolar generalized thermoelastic medium of infinite extent. We are studying the 
propagation of axial symmetric waves which are pure sinusoidal along the axial 
direction. Cylindrical polar co-ordinates (r, 9, z) are considered with z-axis pointing 
upwards. 

We are discussing two-dimensional problem with symmetry about z-axis. 
Therefore, we have 


u = (u r ,0, u s ), 


<j> = (0,<j) e ,0), 



( 1 ) 


where u r , u z and <j> e are related in terms of potential functions (J)' t|/ and f as 


d<|>' d 2 \\i <%' 

w = — + — —,u = — 
r dr drdz 2 dz 



d 

dz 2 


2 \ 


V.4>e =- 


dr 

dr 


( 2 ) 


Making use of eqn. (2) in the field equations derived by Eringen 7 and Green Naghdi 3 
and then eliminating T and T from the resulting expressions, we obtain 
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[□ 2 □ 3 + £ 2 v 2 ] m = o 

and [□, D 4 -T] ^ V 2 ][<)>'] = 0 

dt 


where 


dr 2 r dr dz 2 


.2 


□i=(X + 2p+.K)V 2 - p 


dt 2 ’ 


□ 2 = (p+*)v 2 - p |1 
at 


□ 3 = yV 2 -2^-py^ 


2 9 2 


□ 4 = V^~ T| 


dt 1 


p c* o 2 r 0 

r] = - — , 11 , = 

K * * 


The time harmonic wave solutions of eqn. (3) and (4) are given by 


1 = 


A,KMj) + AK 0 {q,r) 




(3) 

(4) 


(5) 


( 6 ) 


where / = 1,2, ; m = 1, n = 2 and / = 3, 4, ; z?z = 3, re = 4 , 1 stands for \\i, F, (j>' 

T and corresponding to them i = 1, 2, 3, 4 respectively. A, and A\ (z = 1, ....,4) 

are constants and 




(a 2 +a 2 +r) 2 -r\l)+ fcrf -a 2 -t| 2 +r\ 2 2 } +4r) 


, 2_2 

b CT 2 


1/2 




1 rj, 
T1 + — + -T 
. < pc, 


V 


- + - 




4rm, 

pc, 2 


iK' 


(7) 


where 


a 


2 

l 




„ 2 _ 
^2 = 


r 


2K 

Y 
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2 a. 4- 2(1 + K 2 
C \ = 7 5 C 2 
p 



( 8 ) 


The phase velocity of waves in the z-direction is denoted by c (==cd/&) 5 k is the 
wave number and co being the circular frequency K 0 ( ) is the modified Bassel function 
of second kind and zero order. 


A{ and A' (i = 1, ,4), are related as 

A=b_ x A^ A' = b ( A'_ 1} i = 2,4, 

where 


(9) 


b={k l -o\-q])lp, 


pc, 2 

V 

V 

u 


-k' + q) 




( 10 ) 


Boundary Conditions 

The boundary conditions are the vanishing of stresses and temperature on the 
surface of the cylindrical bore i.e. at r = a 


8T n 

t rr -t n = m^= — =0 

dr 


(11) 


where t rn tn, and have been borrowed from the constitutive relations derived by 
Eringen 7 and Naghdi 3 given as 


trr (2|J, d* fC) U rtr — VjT, 


( 12 ) 


trz \^Mr,i d" Wj,r) d" K(tl z r + <j)0), 


(13) 
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m 


re “ 


ftyj , Y^e 

r dr 


(14) 


Using the expression (12) - (14) along with eqn. (2), (6) and (9) in the boundary 
conditions (11), we shall obtain four homogenous equations in four unknows namely 
Au A\, Aj. A non-trivial solution for these unknowns would give the following 
frequency equation 


K t {q t a) nu K 0 {q,a) 

34 A *,M. 

(l^ + K^k 




(a,D, -a,D,)~ 

K,{qa) 


KXq t a) 4 3 K x {q,a) J 


n „ |X(?. fl ) , 1 

— 7 — - + 


Da \!L M,J_] 

K x {q x a) " q x a\ j K, (q,a) T J 


= 0 


(15) 


where 


A =*,<?, 


(P.+y) , .. Ak a ) 

q, a K M, a ) 








,j = 3, 4, 


a, =p(fc 2 + g 2 )+£((>. +# 2 ) ,i = 1,2, (16) 

Eqn. (15) determines the phase velocity of axial symmetric surface waves as a 
functions of dimentionless wave number ka, indicating that the waves are dispersive 
in nature. 

Physically, it must be expected that for every short wave length i.e. for large 
values of ka, the dispersion eqn. (15) should coincide with the dispersion equation for 
the Rayleigh type waves at ta plane boundary of micropolar for the Rayleigh type 
waves at a plane boundary of a micropolar generalized theremoelastic half-space. This 
can be done by putting ka = °o in eqn. (15). For large value of the argument, the 
asymptotic value of modified Bessel functions 12 ' is 
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K 0 (u) = K l (u)' 


I 


■e~“ , («-* oo) 


Hence eqn. (15) reduces to 

{D:b t -D:b 3 ){a l D[-a 1 D , y{2^ + Kyk%b i -qA){D[q l -D[q,) = 0 


where 


z>; = y, d; =(x+ 2^+ *)?; - - vz>, , /= 1, 2; j =3, 4 , 


Special Chases (I) : Neglecting thermal effect i.e. taking v = C = 
and (18), we obtain the following dispersion equations 

%{a x D,-a l D x )-{2\L + K'fk'q\ 



-A^ 


aM 1 1 V 
*.(&«) ? 2 «J. 


= 0 


and 


D;(o 1 o;-o 1 D;)-(2n+x)‘t' ? ,*(A'« 1 -A'?,) -o 


where 


A =(2 ^K)qf 


A~.fag) . 1 

7 ^ / 0 \ 




r 2 \ 

?3 


V 


y 



and 



(17) 


(18) 


(19) 

* = 0 in eqn. (15) 


(20) 


( 21 ) 


( 22 ) 
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Eqn. (20) is the dispersion equation for the propagation of elastic waves in a 
cylindrical bore, through micropolar elastic medium of infinite extent and coincides 
with the eqn. (2.24) of Banerji and Sengupta 10 with the appropriate change in 
notations, Eqn. (21) coincides with the equation of Rayleigh waves at a plane 
boundary of microporal medium obtaining by Kaliski et a/ 16 . 

(II) If we neglect micropolar effect (a = P = y = K-j = 0) in eqns. (15) and (18), we 
obtain the period equation for a generalized thermoelastic medium containing 
cylindrical bore as 




£A' 


go 


B.b. 


^ KM ,«). 


4p 2 Ar A' 2 


KMa) , 1 

Kfaa) q\a 



K x {q t a) 


q A 


-goCiii) 


= 0 


and 

p(F -q':){B[b:-B:b[)-A^kX {q&-q&) = 0 


(23) 


(24) 


where 


g = 2 M ) 


K Xq t a) f 1 
Kjgja) qj a 


+ 




7 = 3,4 


3=[2 Wy 2 +x,fe-* 2 )-v&;], 7 = 3 , 4 , 


i' = 


,'2 


CL 

'a 2 

v j 

* 



= k 2 

1- 


r-* 2 + <?; 


V c 2 7 


7= 3, 4, c" = 


X + 2|i 
P 


♦ 


c 2 



p 


i? 3 _ 4 are again given by equation (7), only with replacement of c\ by c* . 
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(III) Neglecting micropolar and thermal effects together in the medium 
considered, one can obtain the eqn. (2.26) and (2.27) of Banerji and Sengupta 10 and 
represent the dispersion equation obtained by Biot 9 for the propagation of elastic 
waves in borehole. 

Case II-Propagation of waves in a cylindrical bore filled with liquid : Here we 
consider that the cylindrical bore is filled with homogeneous inviscid liquid. For 
axially symmetric waves, the displacement potential <j)o in the liquid satisfies the wave 
equation 


where c\ = XJ p 0 is the velocity of dilatational wave in liquid ; p 0 , are the density 
and bulk modulus of the liquid respectively. The solution of eqn. (26) corresponding 
to the surface wave may be written as 


h = A 5 I 0 {rZ,)e lilb - M) ,foTc<c 3 (27) 

G ) 2 

where C; = k 1 — - and 7 0 ( ) is the modified Bessel function of first kind and order 

zero. Notatingy=c/c 3 , the liquid pressure and radial displacement respectively are (for 
y<l) 

P = Po® 2 4 4 [r*(l - / T 1 e '^' ] (28) 

q r = k(l - y ! )' /2 A S I, [r*(l - / f ] e‘ ( ^ u) (29) 


Boundary conditions : The appropriate boundary conditions at the interferace r = a 
between liquid and micropolar generalized thermoelastic medium are 
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The first two boundary conditions can be written as 

- P _ L 

I,-'.' 

which is a composite boundary condition to the matching of mechanical impedance. 
The last three boundary conditions are same as in Case I. These three boundary 
conditions and boundary condition (3 1) with the use of eqns. (6), (9), (12) - (14), (28), 
(29) give four equations in four unknowns namely A u A\, A 3 and A\. Fro non-trivial 
solution of these four homogenous equations, we must have the determinant of the co- 
efficient of the unknowns equal to zero, which yields (fory <1). 



^o(<?4a) 

K,{q,a) 


D 4 b 3 


^p(?3 a ) 

*>M 


(a,D 2 -<7 2 A)~ 


(2 ix + Kfk 2 



K oia^) „ h ^o(?3 a ) 
3 


D s =2 


(32) 


where 


D] = £>, - mq j j = 3, 4, 


D 5 = 





MD X 

( 2\i+K ) 


M = 


P 0 C0 2 I 0 


*(i -y'T /,! 

ko-yf] 


(33) 


If thermal effect is neglect then it can be verified that dispersion eqn. (32) recuces to 
that obtained by Banerje and Sengupta 10 for the relevant problem in micropolar elastic 
medium. 
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In the absence of microrotation effect , eqn. (32) reduces to the dispersion 
equation for surface wave propagation in a generalized thermoelastic medium 
containing a cylindrical bore hole (filled with inviscid liquid) as 






b'A 


Kq (q 3 g) 


4[i 2 k 2 q[ 2 


K\{q[a) q[a 


+ — Mq[ 


q]b: 


, K 0 {q 4 a ) 



(34) 


where 


Bj = Bj - Mq } 7 = 3,4, 

Numerical Results and Discussions 

In order to study the problem numerically, frequency eqn. (15) and (32) are solved for 
phase velocity C 0 = (c/ci) by giving different values of wave number in non- 
dimensional form ( ka ) to study the thermal and microrotation effects on dispersion 
equation, we take the case of magnesium crystal 17 like material subject to thermal 
disturbance for numerical calculations. The physical constants used are 

p = 1.74 g/cm 3 ,7 = 0.2 x 10 ~ 15 cm 2 , X = 9.4 x 10 n dyne/cm 2 

p = 4.0 x 10 11 dyne/cm 2 , K = 1.0 x 10 n dyne/cm 2 

y = 0.779 x vr* dyne, K* = 1.0005 cal dyne/g °C cm 2 

C* = 0.23 cal/g °C, e = 0.073, v = 0.89 dyne/cm 2 °C, 


T 0 = 23 °C. 


The following values are taken from Ewing et a/. 18 for inviscid liquid : 
Ao = .214 x 10 u dyne , p 0 = 1.0 g/cm 3 . 
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Fig. 1 -variations of phase velocity with wave number. 
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For the above values of relevant parameters, we have calculated the normalized 
phase velocity C 0 against the values of normalized wave number ka, as depicted in 
Fig. 1 . It is found that there exist infinite number of modes of propagation of surface 
waves for various values of wave numbers; the variations of only fundamental modes 
are presented graphically. These variations of phase velocity in micropolar 
generalized thermoelastic (MGTE) medium, generalized thermoelastic (GTE) medium 

and micropolar elastic (ME) medium are shown by (a) solid line ( ), small dashed 

line ( ) and long dashed line ( ) respectively in case of empty bore, (b) solid 

line with centre symbols (— x— x— x— ), small dashed line with centre symbols (-x -x -x) 
and long dashed line with centre symbols (-x -x -x Respectively in case when bore is 
filled with inviscid liquid. 

From Fig. 1, it is observed that in MGTE medium, the phase velocity decreases 
gradually with the increase in wave number as depicted by the two curves, one solid 
and other with centre symbols, also indicating that the effect of liquid is very small in 
this medium , If, micropolarity is neglected i.e. in GTE medium, the values of phase 
velocity are large but the variations lie in a very small range in comparison to MGTE 
medium as can be seen by the two curves; one small dashed and other with centre 
symbols, also characterising the significant effect of liquid in GTE medium. In the 
absence of thermal effect i.e. in ME medium, a similar trend i.e. decrease in phase 
velocity with the increase in wave number is observed and the values of phase 
velocity are large in comparison to MGTE medium. Since the difference in the values 
of phase velocity in ME medium for the cases of empty bore and liquid filled bore is 
very minute in the range 3.6 < ka < 9.6, therefore the effect of liquid is not visible 
from Fig. 1 in this range. 

One of the interesting phenomenon observed on comparison of the various curves 
in fig. 1 is that the values of phase velocity become large after the removal of 
microrotational or thermal effects from the medium, although the behaviour of 
variations is similar in all the curves. Also, the effect of inviscid liquid present in the 
bore hole on propagation of surface waves is more significant in ME and GTE 
medium in comparison to MGTE medium. This effect can be seen on comparing the 
three curves with centre symbols and the corresponding three curves without centre 
symbols in Fig. 1. 
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Abstract 

In this paper we study a special type of topology and its property in the plane, called Sun 
Topology. 

(Keywords : sun topology/density topology) 

Introduction 

We know that there are three well known topologies in the plane (Euclidean 
topology, Fine topology, Usual density topology) are linearly ordered. The fine 
topology, is strictly finer than Euclidean topology, and has several useful geometric 
properties which disappear when use switch to the finer density topology. 

Definition : The Lebesgue measurable set A in the complex plane is said to be sun 
open if every x e A for almost every a e [0, 2n] 3 ? a > 0 : {z e C : z = x + 1 (cosa + i 
sina), |r| <t a } <^A-. 

Lemma : The density topology is strictly finer than the sun topology. 

Proof : Let a set A be sun open at x. For each r > 0 denoted by B r the set of all p e 
[0, 2 tc] for which {zeC\ z-x + t (cosp + i sinP), |t| <r} c A. We see that X(5 r ) — > 
2n for 0, hence 
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Hence x is a point of density of A. The set A density open at x. Hence the density 
topology is finer than the sun topology. Again let the set C — {z e C : l/|z| e N} is 
density open at 0 and is not sun open at 0. Hence the density topology is strictly finer 
than the sun topology. 

Result : Now we introduce the fine topology from potential theory in the plane. If 
we denote for an arbitrary point* e C, a set A c C and n e iVthen 

A„ = {z e C\A, 1/2” < |z-x|, 2/2”}. 

We can characterize points * at which a set is finely open as those for which the 
series 


x — V) 

y 

log Cap [t„(x)] 


converges (by Wierer’s test), where Cap denotes the outer logarithmic Capacity [1, 
see theorem 9, 10]. 

Theorem : The sun topology is strictly finer than the fine topology. 

Proof : Let us consider 0 e A and A is finely open at 0. We now restrict to K = {z 
= x + iy e C:*>0,y>0} and we show that for almost every a e [0, tt/2] a t a > 0 : {z 
e K:z = t (cosa + i sina), 0 < t < t a } c A. for n e N. 

let 


A„ = {z e K-A, 1/2" < \z\ < 2/2"}. 

Let B„ =p(A„), where p is the projection from the origin on the line through the points 
[1/2”, 0] and [0, 1/2"] also let C n = R(A-n), where R is the projection from the origin on 
the circle {z e C : |z| e 1/2"}, then the outer logarithmic capacity increases under 
contractions, now Cap(5„) < Cap (A„) for the measurable liner set B„ we can use the 
estimate. The linear Lebague measure on the line through the points [1/2", 0] and [0, 
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1/2"] by X(B n ) < 4 Cap(J3„) [2, see page 173], There is a simple Lipschitz mapping 
from B„ -» C„. Hence the linear Lebesgue measure of C„ on the circle {z e C : \z\ = 
1/2"}, is estimated by A ,(C n ) < 2\(B n ). The set is finely open at 0, hence the series 


V 

& logCap(^„) 


converges. 


Now we conclude that the series 


00 — Yl 

SlogA(C„) 


converges. 


We know that A(C„) e [0, 7t/2" +1 ]. Each point at C N stops the radius from the 
origin at the distance 1/2". The measure of the angles which stop at C„ is q„ - 2"A(C„). 
We show that converges. We see that the measure of angles in [0, 2/n] which 

stop at same distance due to C„ (Also A„} for n > m is smaller thanE” =1 ^ H . Hence 
arbitrarily small for m big enough. Hence almost every angle at [0, 2 /tc] stops at same 
positive distance. The set A is sun open at 0. The set C-R + is obviously sun open at 
origin and is not finely open at the origin due to the fact that each finely open set 
contains with any point arbitrarily small circles centered at the point [3, see theorem 
10.14], The sun topology is strictly finer than the fine topology. 

Remark : The sun topology under the name the core topology, almost everywhere, 
was already studied byHorbaczewska [4, see page 416], it has quite recently been 
proved by Wagner-Bojakowsak and Wilczynski [5, see theorem 2, page 363]. 
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Abstract 

A tetrahedral finite element based mathematical model has been developed to obtain the effects 
of anisotropy on the temperature distribution (TD) in three dimensional regions of human skin 
and subcutaneous tissues (SST). Galerkin's finite element technique has been used to compute 
steady state nodal temperatures for a comparative study of the effects of anisotropy, isotropy, 
blood flow, metabolic heat generation and evaporation by considering eleven different 
combinations of these parameters. It has been found that fall of temperature due to evaporation 
is more in isotropic medium than that of in an isotropic medium. The study also attributes that 
not only evaporation brings down the temperature of SST region, but similar role is also played 
by the components of thermal conductivity taken highest in the direction of width or length of 
SST. The study may be useful in formulating sophisticated mathematical models for the heat 
flow problems with tumours. 

(Keywords : anisotropic skin/blood flow/tetrahedral element/Galerkin’s finite element.) 

Introduction 

All animals/human beings live in variable thermal environment with varying 
internal metabolic heat generation. They are able to maintain a stable internal 
temperature at level which is generally above that of the environment. The thermal 
gains and losses vary with different parameters participating at different levels and 
times. Heat migration in human skin and subcutaneous tissues (SST) is the 
consequence of blood flow, metabolism and environmental conditions. Heat 
generation and its loss to the environment are the continuous processes in human 
body. Blood flow plays a vital role in carrying excessive heat from the body to the 
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skin surface where it is given off to keep the core temperature constant within narrow 
limits. The skin consists of two layers known as epidermis and dermis. The outer dead 
layer of epidermis, the stratum comeum, is composed of flattened cells. The dermis, 
which is the deeper and thicker of the two major layers of the skin, is composed of 
matted masses of connective tissues and elastic fibres through which pass numerous 
blood capillaries and nerves. The subcutaneous tissues, which attach the skin with the 
rest of the body, contain fat cells, connective tissues, larger blood vessels, lymph 
vessels and nerves. The thermal characteristic of tissues plays a very important role in 
heat regulation in SST subjected to external surrounding conditions and variable 
metabolic activity. 

Various mathematical models for computing temperature distribution in SST 
regions have been developed by many investigators. Pennes 1 formulated an energy 
equation for living human tissues by idealizing the human forearm as a cylinder 
excluding the axial and angular effects. Gordon and Roemer 2 studied the effects of 
radial node spacing on finite difference calculations of temperature in living tissues. 
The thermal recovery of skin temperature after cooling has been studied by Stekette 
and Vender Hoek 3 . The SST was considered as a single region in these models with 
constant values of biophysical parameters like blood flow, metabolic heat generation 
and thermal conductivity. These models were oversimplified version of the 
multilayered nature of SST region with variable biophysical parameters. Saxena 4 , 
Saxena and Pal 5 , and Pal and Pal 6 ' 8 developed improved multilayered one dimensional 
(1-D) mathematical models with heterogeneous variation of biophysical parameters 
and studied the temperature pattern in human SST under various environmental and 
physiological conditions. 2-D mathematical models for the steady state temperature 
distribution of SST region with triangular finite elements were developed by Arya and 
Saxena 9 and Pal and Pal !0 . The two dimensional mathematical model 10 revealed that 
the temperature profiles reflect the dependence of temperature distribution not only on 
the environmental conditions and biophysical parameters but also on the geometry of 
the SST region. Pardasani and Adlakha 11 investigated steady state temperature 
distribution in annular tissue layers of a human or animal body by applying cubic 
spline technique. In a later work, Pardasani and Adlakha 12 used coaxial elements to 
study radial and angular heat distribution problem in human limbs by considering core 
temperture variable in angular direction. It has been found that coaxial elements 
reduce computational work of the finite element technique in comparison to triangular 
elements. All these studies are confined to one and two dimensional models only. The 
difference between various multilayered mathematical models was of the number of 
layers and sublayers into which SST was divided according to the geometry and 
behaviour of biophysical parameters. Crisp 13 and Stolwijk 14 have discussed the 
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properties of SST region and suggested that it exhibits anisotropic behaviour. This 
aspect has not been addressed by the previous investigators. To understand the 
complete effect of anisotropy on temperature distribution, three dimensional 
mathematical models of SST region are required to be formulated. 

In the present work, a 3-D mathematical model of anisotropic SST region has 
been developed by dividing the region into tetrahedral finite elements. Heterogeneous 
variation in biophysical parameters like blood flow, metabolic heat generation and 
thermal conductivity have been considered for computing the steady state nodal 
temperatures by Galerkin's finite element technique. Eleven different combinations of 
anisotropy, blood flow, metabolic heat generation and evaporation have been 
considered to study and compare their effects on temperature distribution in three 
dimensional SST. Such studies may be useful to medical scientists in formulating 
advanced mathematical models for the detection and diagnosis of tumours. 

Mathematical Model 

It T denotes the temperature at any point (x,y,z) of SST region which is assumed 
to be anisotropic, the governing bioheat equation for the steady state case is 

k x (d 2 T/^ 2 )+k y (d 2 T/d y 2 )+k z (d 2 T/a: 2 )+w b Cb(T b -T)+Q =o (i) 

where k Xt k y% k z = thermal conductivities of the tissues in x(width), y(length) and z 
(thickness) directions respectively; w b = rate of blood mass flow; c b = specific heat of 
blood; Q = rate of metabolic heat generation and 7* = temperature of blood when it 
enters the tissue compartment. Let a and b the width and length respectively of SST 
under consideration. The thicknesses of epidermis, dermis and subcutaneous tissues 
are considered to be C\ C 2 , c 3 respectively. The total thickness 'c' of SST can be 
expressed as 

c=ic, p) 

/'-l 

If R, D and S,(i=\-6) are the region, domain and the boundary respectively, then 

R=D+ts, 

1=1 


( 3 ) 
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where S'/ : z = 0; S 2 : z = c, S 3 : y = 0, S 4 : y = b£p -Q-,S(,:x = a. 

Each layer of SST has been divided into six, eight cornered, blocks by trisecting 
and bisecting its length and width respectively (Fig. 1). Each block is further divided 
into five tetrahedral elements’ 5 as shown in Fig. 2. Thus, whole of SST has been 
discretized into 90 tetrahedral elements which are connected to each other by forty 
eight nodes. Details of these blocks and their respective node and tetrahedral element 
numbers have been given in Table 1 . 

Table 1- Blocks, their respective node numbers and tetrahedral element numbers 

Block Nos. Node Nos. Tetrahedral Element Nos. 


1 

1 

2 

5 

4 

7 

8 

11 

10 

1-5 

2 

4 

5 

14 

13 

10 

11 

20 

19 

6-10 

3 

13 

14 

17 

16 

19 

20 

23 

22 

11-15 

4 

2 

3 

6 

5 

8 

9 

12 

11 

16-20 

5 

5 

6 

15 

14 

11 

12 

21 

20 

21-25 

6 

14 

15 

18 

17 

20 

21 

24 

23 

26-30 

7 

7 

8 

11 

10 

25 

26 

29 

28 

31-35 

8 

10 

11 

20 

19 

28 

29 

38 

37 

36-40 

9 

19 

20 

23 

22 

37 

38 

41 

40 

41-45 

10 

8 

9 

12 

11 

26 

27 

30 

29 

46-50 

11 

11 

12 

21 

20 

29 

30 

39 

38 

51-55 

12 

20 

21 

24 

23 

38 

39 

42 

41 

56-60 

13 

25 

26 

29 

28 

31 

32 

35 

34 

61-65 

14 

28 

29 

38 

37 

34 

35 

44 

43 

66-70 

15 

37 

38 

41 

40 

43 

44 

47 

46 

71-75 

16 

26 

27 

'30 

29 

32 

33 

36 

35 

76-80 

17 

29 

30 

39 

38 

35 

36 

45 

44 

81-85 

18 

38 

39 

42 

41 

44 

45 

48 

47 

86-90 
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Numbers indicate node numbers. 
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It is assumed that (i) Si coincides with the body core (ii) S 2 is exposed to the air 
environment and loss of heat is due to convection, radiation and evaporation of sweat 
(iii) surfaces Si (i = 3-6) are insulated. Thus, we have 


T(x, y, 0) = 37 °C, 

(4) 

k z (dT/dz)=h(T-T a ) + E, at z = c, 

(5) 

where h = heat transfer coefficient, E = evaporative heat loss, 


h y (dT/dy) = 0 on S 3 and S 4 , 

(6) 

k x {dT/dx) = 0 on S 5 and S 6 , 


The evaporative heat loss E is calculated from the following equation 


E = X. rate of evaporation of water, 

(7) 

where X = latent heat of vaporization of water. 

The trial temperature 7^ within a typical tetrahedral element 'e\ whose node 
identifiers p(x p y Pt z p ), q(x q _ y q z q ), r(x r y r z r ) and y Si z s ), are taken in a 

counterclockwise direction, is assumed to be 

f e) = N (e) (x ,y,z) T (e) , 

( 8 ) 

where 


N (e) (x,y,z) = [N p (x,y,z) N q (x,y,z) N r (x,y,z) N s (x,y,z)], 

(9) 

T (e) = [ T p T q T r T s ] 

( 10 ) 


and T (e) are shape function and nodal matrix respectively, Prime donotes the 
transpose of the vector. 

Substitution of trial solution ( 8 ) into equation ( 1 ) (after rewriting it for the element 
'e') results into a residual. The minimization of this residual with the help of 
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orthogonality conditions with respect to the basis function produces the following set 
of integrals. : 

I (t) ={k[ e) {8 2 T (e) /8x 2 ) + k ( y e) (d 2 T (e) / 8y 2 ) + k ( * ) (8 2 T {e) /& (2) ) 

+ w ( b e) c ( b e) (T b -T^) + Q (e) ,N j )= 0 (11) 

where j = p, q, r, s. Green's theorem converts the equation (11) into the following 
form : 

- (k'f (8T (e) / etc), (8Nj / dx)) - (k ( f> {8T ie) / 8y), ( 8Nj / 8y)) - 

(k[ e \8T (e) /dz),(dN J /dz))-(w ie) T {e \N J )+ J |[^ e) (5r (e) /8x) + 

& 

k<f ] (8T (e) / 8y] + k ( ; ] ( 8T (e) / 8z)]N jds (e) + (w (e) T b ,Nj) + (g (e) , Nj ) = 0 

(12) 

where w (e) = w { b e) c ( b e) 

Combining (8), (9) and (12) and imposing the Cauchy boundary conditions (5) 
and (6), we get 

-(fcfT (e) (aJV,. IdxUdNj /dx))-(k^\dN 0] /dy),(dN J /dy))- 
< T (e) )5A r i /&), (SAL / 5z)) + (w (e) r A (w (e) T (e) AL , AL ) + (2 (e) , AL ) - 

J JN - T a ) + E]ds^ =0; 

,00 

i,j=p,q,r 9 s. (13) 

The above equation can be written in the matrix form as 
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i K, F<« = £pf , (14) 

m=l m=l 

where 

[K[ e) ] = (k[ e) (8N I / dx), (dN J /dx)) + (kf (dN, / dy), (« dN] / dy)} + 
(k { / ) (dN i / dz),(6Nj / &)}; 

[k 2^ ] = (w (e) N^Nj ); [K?]= j <fc< e) ;P, <e) =(w< e) T b ,Nj)-, 

4 e) 

P 2 (e) = (2 (e) , TV,); P 3 (e) = J -£]#,&<*>; /, j=p, q, r, s. (15) 

,(«) 

The function Nj is unity at the no 'j' and zero at all the other nodes and can be 
expressed as 

Nj(x,y,z)=(aj+bjX+Cjy+djZ)/6V < - e) , j=p,q,r,s, (16) 

where 

1 */> z p 

y(e) = x q yq z q =v 0 i ume 0 f tetrahedron pqrs, (17) 

1 x,. y r z r 

1 y s z s 

X q yq Z q 1 ^ Z q ' 

a p = x r y r z r , b p - - 1 y r z r 

x s y s z s 1 y s z s 

* q \ z q x q y q 1 

C p =-x r 1 z r , dp=-x r y r 1 (18) 

1 *, y s i 
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Other constants are obtained through a cyclic permutation of subscripts p,q,r and s. 
Equations (17)-(18) are valid only when the nodes p,q and r are numbered in a 
counterclockwise manner when viewed from the node 's'. The simplification of the 
above integrals yields the following matrices for the element 'e' : 


[K' e) ]=(l/36F (e) ) 


8n Sn 8n 8u 

Sl2 8 23 8 24 


symmetry 


•§33 8 34 


8 44 


[Kf]=(w/ e) F (e) /12) 


2 1 
1 2 
1 1 
1 1 


1 1 
1 1 
2 1 
1 2 



1 

1 

2 

0 


0 

0 

0 

0 


P, (e) =(w^ e) T b V (e) / 4)[1 1 1 1]'; 


P^ e) = (0 (e) F w / 4)[l 111]'; P 3 (c) = [(hT a - E)A%\ /3][l 1 1 0]' ; 


gn=rt%+K?c 2 p+ k?d 2 p y> 

8i2=^%b q+ k^c p c q+ k^d p d ti y, 

gn=&Xb r+ k?c p c r+ k?d p d r y 
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gn = (*X + + *X); &, = (tfty, + *? V, + 

ft, = + k T c , c :. + k ? d , d 'X ftj = <k[''b: + k’’‘c‘ + ); 


ft, =(*X, +*JV,+*X</.); ft. =<*X +*X +*X); 


(19) 

and is the surface area of the face pqr. Each element contribution is computed 

by using element equation (14) and assembly of these contributions for all the ninety 
elements of the SST region gives the following system matrix equation : 

90 3 r 1 r 1 90 3 

Il[KL e) ][T (e) ]=ZZPi e) (20) 

e~\ m = 1 e=l m-\ 


Above equation is solved after the insertion of Dirichlet condition (4), to get the nodal 
temperatures. 


Results, Discussion and Conclusions 

. « 

The chief interest of the present investigation is the comparative study of the 
effects of anisotropy, rates of blood flow and evaporation on temperature distribution 
in three dimensional regions of human SST. The width, length and thickness of SST 
can be assigned any value depending on the sample under study. In the present case 
we have considered a= 1.00 cm, b- 0.9 cm, c= 0.7 cm, ci=0.15 cm, C2=0.30 cm, C3=0.25 
cm. Eleven different combinations of biophysical parameters (elementwise) and 
evaporation denoted by C/(z=l-l 1), given in Table 2, have been considered to compute 
the nodal temperatures at environmental temperature T a ~ 25 °C. Combinations 
C,(/=l,2,3) are considered for isotropic medium. Evaporation rate is considered higher 
in C 2 than that of C\. The rates of blood flow and metabolic heat generation have been 
considered to be double in C 3 than that of C\ (or C 2 ). Combinations C,(z=4-1 1 ) have 
been taken for anisotropic medium. C 4 , C 5 , C 6 incorporate anisotropy in epidermis, 
dermis and subcutaneous tissues respectively whereas C, (r=7-l 1) consider anisotropy 
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Table 3- Computed nodal temperatures in three dimensional regions of human SST. 


Node 

Nos. 

c, 

c 2 

C 3 

c 4 

C 5 

C 6 

Cl 

C 8 

c 9 

Cio 

C u 

1 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

2 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

3 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

4 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

5 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

6 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

7 

35.13 

33.35 

35.30 

35.09 

34.97 

35.74 

35.59 

33.85 

34.25 

34.85 

35.68 

8 

35.66 

34.37 

35.81 

35.64 

35.63 

36.22 

36.18 

35.16 

34.74 

35.32 

36.26 

9 

36.11 

35.21 

36.23 

36.09 

36.13 

36.60 

36.62 

36.10 

35.18 

35.71 

36.67 

10 

35.57 

34.19 

35.72 

35.54 

35.50 

36.14 

36.07 

34.92 

34.67 

35.21 

36.15 

11 

35.66 

34.37 

35.80 

35.64 

35.62 

36.20 

36.16 

35.11 

34.74 

35.33 

36.23 

12 

35.77 

34.57 

35.91 

35.75 

35.75 

36.30 

36.27 

35.33 

34.82 

35.44 

36.34 

13 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

14 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

15 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

16 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00* 

37.00 

37.00 

17 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

18 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

37.00 

19 

35.78 

34.58 

35.91 

35.75 

35.76 

36.30 

36.28 

35.36 

34.83 

35.44 

36.35 

20 

35.67 

34.39 

35.81 

35.65 

35.63 

36.20 

36.17 

35.15 

34.74 

35.34 

36.24 

21 

35.57 

34.20 

35.72 

35.54 

35.50 

36.15 

36.08 

34.93 

34.67 

35.21 

36.15 

22 

36.11 

35.22 

36.23 

36.09 

36.14 

36.60 

36.62 

36.10 

35.18 

35.71 

36.68 

23 

35.67 

34.40 

35.81 

35.65 

35.64 

36.22 

36.19 

35.20 

34.75 

35.32 

36.27 
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Table 3 Contd.. 


24 

35.14 

33.38 

35.31 

35.10 

34.98 

35.75 

35.59 

33.88 

34.26 

34.85 

35.69 

25 

32.45 

28.23 

32.66 

32.35 

33.27 

32.98 

33.78 

29.94 

29.82 

31.56 

33.89 

26 

33.32 

29.86 

33.54 

33.26 

34.22 

33.84 

34.72 

31.93 

30.65 

32.42 

34.83 

27 

34.03 

31.16 

34.25 

33.96 

34.98 

34.53 

35.45 

33.46 

31.38 

32.97 

35.56 

28 

33.11 

29.47 

33.33 

33.02 

33.99 

33.64 

34.49 

31.44 

30.54 

32.10 

34.60 

29 

33.37 

29.96 

33.59 

33.32 

34.25 

33.88 

34.76 

32.03 

30.72 

32.45 

34.86 

30 

33.46 

30.09 

33.68 

33.40 

34.36 

33.97 

34.86 

32.18 

30.81 

32.48 

34.97 

31 

30.84 

25.23 

31.04 

31.39 

31.63 

31.34 

32.79 

27.85 

26.92 

29.47 

32.90 

32 

31.46 

26.30 

31.67 

32.14 

32.30 

31.95 

33.55 

29.39 

27.36 

30.10 

33.66 

33 

31.92 

27.05 

32.13 

32.67 

32.81 

32.40 

34.12 

30.48 

27.75 

30.34 

34.23 

34 

31.23 

25.92 

31.44 

31.88 

32.06 

31.74 

33.30 

28.88 

27.27 

29.72 

33.41 

35 

31.54 

26.45 

31.75 

32.22 

32.38 

32.03 

33.62 

29.54 

27.48 

30.14 

33.73 

36 

31.58 

26.45 

31.79 

32.28 

32.44 

32.07 

33.70 

29.60 

27.52 

30.07 

33.81 

37 

33.48 

30.14 

33.69 

33.42 

34.38 

33.99 

34.88 

32.26 

30.83 

32.48 

34.99 

38 

33.40 

30.08 

33.62 

33.36 

34.28 

33.91 

34.79 

32.18 

30.75 

32.47 

34.90 

39 

33.13 

29.52 

33.35 

33.04 

34.00 

33.65 

34.50 

31.50 

30.56 

32.11 

34.62 

40 

34.03 

31.18 

34.25 

33.97 

34.98 

34.53 

35.45 

33.45 

31.40 

32.97 

35.56 

41 

33.35 

29.97 

33.57 

33.29 

34.24 

33.87 

34.75 

32.04 

30.68 

32.45 

34.86 

42 

32.47 

28.28 

32.68 

32.37 

33.29 

32.99 

33.79 

29.97 

29.85 

31.57 

33.90 

43 

31.61 

26.54 

31.82 

32.31 

32.47 

32.10 

33.74 

29.73 

27.56 

30.07 

33.84 

44 

31.66 

26.92 

31.87 

32.32 

32.50 

32.15 

33.73 

29.98 

27.63 

30.25 

33.84 

45 

31.25 

25.96 

31.46 

31.90 

32.07 

31.75 

33.31 

28.92 

27.31 

29.73 

33.42 

46 

31.92 

27.01 

32.13 

32.67 

32.80 

32.40 

34.11 

30.41 

21.71 

30.34 

34.22 

47 

31.50 

26.43 

31.71 

32.17 

32.34 

31.99 

33.59 

29.51 

27.41 

30.15 

33.69 

48 

30.85 

25.27 

31.06 

31.41 

31.65 

31.36 

32.80 

27.88 

26.95 

29.48 

32.91 
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through out the region i.e. epidermis, dermis and subcutaneous tissues. The rate of 
evaporation in C 7 and C% have been taken equal to those of C , and Ci respectively. In 
C 7 the component of thermal conductivity (CTC) is considered highest in z direction 
than that of x and y directions. C9 considers CTC highest in x direction than that of y 
and z directions whereas C10 has CTC highest in y direction than that of x and z 
directions. In C\\, the rates of blood flow and metabolic heat generation have been 
assumed to be double of those of C 7 . All the nodal temperatures computed for these 
eleven combinations have been listed tin Table 3 . 

Comparative study of cases Ci with C 7 (isotropic medium) and C 7 with C% (fully 
anisotropic medium) reveals that higher rate of evaporation results in fall of the 
corresponding nodal temperatures (CNT) in both the media. It has been observed that 
the effect of evaporation is maximum at the skin surface and this effect decreases 
towards the body core {i.e. decreases with respect to depth). The difference between 
the CNT of Ci, C2 is more than that of C 7 , C 8 throughout the region. This depicts that 
the fall of temperature due to evaporation is more in isotropic medium than that of in 
anisotropic medium. 

Comparisons of results of C, with C3 and C 7 with C\\ show that there is an 
increase in CNT due to increase in blood flow and metabolic heat generation in 
isotropic as well as in anisotropic media. The study implies that the heat carried by 
blood flow to the skin surface is reduced by the anisotropy of the region. The results 
also infers that the effect of evaporation on TD is more dominant than that of blood 
flow and metbolic heat generation in both types of media. 

Comparisons of the results of C 7 , C9, C10 (fully anisotropic) with Ci (isotropic) 
reveal that the CTC, taken highest in z direction (C 7 ) elevates the temperature 
throughout the SST region whereas the same taken highest in x direction (C 9 ) or in y 
direction (C 10) results in the fall of temperatures from that of C\. The decrease in 
temperature is more in C9 than that of in C10. In all these cases the effect is maximum 
at the skin surface and their respective effect decreases with depth. It implies that 
there is a correlation between the variation of temperature and the CTC in the 
direction of length, width and depth of SST. This attributes that not only evaporation 
brings down the temperature of SST region, but similar role is also played by CTC 
when it is taken highest in the direction of width or length of SST. 

Comparison of nodal temperatures for C,(/= 4 - 6 ) with C\ reveals that the 
incorporation of anisotropy in epidermis, dermis and subcutaneous regions separately 
increases the skin surface temperatures in general, but the increases is highest in C5, 
higher in C 4 and high in C 6 than that of C\ (isotropic case). It is interesting to note that 
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the anisotropy of epidermis increases the skin surface temperatures only whereas the 
anisotropy of dermis increases the temperature of skin surface and of the interface of 
epidermis — dermis. The temperature of whole of SST region is elevated due to the 
anisotropy of subcutaneous tissues. However, their respective effects are nearly 
uniform at skin surface, interface between epidermis-dermis and interface between 
dermis-subcutaneous tissues. It is important to note that in all the three cases (C4, C5, 
C 6 ) the CTC has been taken highest in the direction of depth. These results are of 
interest and may be useful in the formulation of three dimensional sophisticated 
mathematical models for the heat flow problems with tumours in animal/human body. 
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Abstract 

Let (. X,d) be a metric space and T : X -+ X a self-mapping of X. In metric fixed point theory 
there are many fixed point theorems for contractive-type or expansive-type mappings. In this 
paper we define a new class of mappings which besides contractive mappings contains some 
expansive mappings. This class is defined by the following condition: 

min {hd(Tx,Ty\ (d(x, Tx)+d(y, Ty))} < (2 k+a) max {d(x,y), — m(x, y)} 

k 

where k > 0 is arbitrary and 0 < a < 2, h > 2k +a and m (x,y) = min {d{x,Ty), d(y,Tx)}. In 
complete metrically convex metric spaces these mappings have a fixed point. 

(Keywords : metrically convex space/fixed point). 


Introduction 

Let ( X,d) be a metric space and T a self-mapping on X. A solution of the equation 
Tx = x is called a fixed point of T. A question when this equation has a solution 
generates a theory, which began in 1912 with the work of Brouwer 1 , who proved that 
any continuous mapping T of an «-ball into itself has a fixed point, and was followed 
in 1922 by Banach’s Contraction Principle, already obtained in particular situation by 
Liouville, Picard and Goursat, which states that any mapping T of a com-plete metric 
space X into itself satisfying the inequality 

d(Tx,Ty ) < Xd(x,y), 0<X<\ (1) 

for all x,y in X, has a unique fixed point. Schauder 2 , Browder 3 and many others have 
added to and generalized these basic results for continuous mappings. 
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In 1973 Wong 4 proved some fixed point theorems for operators, not necessarily 
continuous. A class of contractive type mappings, not necessarily continuous, were 
significantly extended and generalized for single-valued mappings (5-15), as well as 
for multi-valued mappings (16-18). 

One of the most general contractive-type definitions for which fixed point 
theorems have been proved, as pointed out by Rhoades 13 , is introduced in (8) and 
states: 

d(Tx,Ty ) < Xmax {d(x,y), d(x,Tx), d(y,Ty), d(x,Ty), d(y,Tx)}, 0<A.<1 (2) 

However, in many applications, the mapping involved is not of contractive type 
and does not satisfy anyone of the contractive type definitions. 

The purpose of this short paper is to define and to study a class of mappings, not 
necessarily continuous, which on some parts of domain may be of contractive type, 
and on other parts may be of expansive type. For that class we shall prove fixed point 
and common fixed point theorems. One example is presented to show that this class 
contains mappings which are not necessarily contractive, nor expansive on a whole 
domain. 


Main results 

We need the following definition of metrically convex metric space due to 
Menger. 

Definition 1 : (19, p.41) A metric space (X,d) is said to be metrically convex if for 
each x,y in X, x*y, there is a z * x, z * y for which 

d(x,y) = d(x,z ) +d(z,y). 

A subset CcTis said to be metrically convex if subspace (C,d) is metrically 
convex. 

Lemma 1 : (Menger) : If X is a complete metrically convex metric space, then for 
any A.,0 < X < 1 and any x,y in X, there exists z in X such that 

d(x,z) = Xd(x,y) and d(z,y) = (\-X) d(x,y) 


Proof: See (19, p.41, Theorem 14.1). 
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Now we shall prove the following result. 

Theorem 1 : Let C be a nonempty closed and metrically convex subset of a 
complete metric space X and let T : C -> C be a self-mapping on C satisfying 

min {hd(Tx,Ty), ( d(x,Tx ) + d(y,Ty))} <(2 k+a) max {d(x,y), — m(x,y)} (3) 

k 

for all x, y in C, where k > 0 is arbitrary, a and h are such that 

0 < a < 2, h > 2k +a (4) 

and 

m(x,y), = min {d(x,Ty), d d(y,Tx)}. 

Then Thas at least one fixed point in C. 

Proof : Let Xo in C be arbitrary. Since C is convex, by Lemma 1 we may choose X\ in 
C such that 

d{x 0 , ) = -^~d{x Q , Tx 0 ) and d{x v Tx 0 ) = ~^—d(x Q , Tx 0 ) , 
k+\ k + 1 

or rewritten, 

d{x () , Tx 0 ) = (k + \)d(x 0 , x,), d(x p Tx 0 )=kd(x 0> x,). 

Again by Lemma 1 we may choose x 2 in X such that 

d(x l ,TX])-{k + \)d(x v x 2 ), d(x 2 , Tx l )=kd(x v x 2 ). 


Continuing in this fashion we can choose the sequence {x„} in X such that 
d{x „_ „ Tx nA ) = {k + \)d{x nA , x„), 


( 5 ) 
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d(x n> Tx n _\ ) - kd{x n _ j, x n ) . 

We shall show that {x n } is a Cauchy sequence. Setting x = x„-i, y = x n in the 
definition inequality (3) yields 

mm{hd(Tx n _ v Tx„), [d[x „_ ,, Tx nA )+d{x„, Tx „ )]} 




Using (5) we get 

min {hd(Tx n _ v Tx n ), (k + 1) [d{x nA , x n ) + d(x„, x„ + ,)]} 



Suppose at first that from (6) we have 

hd(Tx„, v Tx„)<{2k + a)d(x nA , x„). (7) 

By the triangle inequality we have 

| 14 * 1 ,. Tx „A- d ( Tx n-\’ Tx n)- 

Hence, using (5), we get 

\{k + \)d{x n , x n+l )-kd{x n _ v x n \<d{Tx nA , Tx „). (8) 

If (A: + 1) d{x n , x, ;+1 ) > kd x „ ) , then from (7) and (8) 

h(k + \)d{x n , x n+l )-hkd(x n _ x n )<(2 k + a)d{x nA , x„) 
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and hence 


,/ \ hk + 2k + a .( \ 

d\x n , x n+l ) < d(x n _ v x„ ) . 


( 9 ) 


If {k + 1 )d(x n , x„ +t ) < kd(x„_ v x n ), then this immediately implies 


d{x n , x n+x )<- — -d{x n _ v x„). 

K+ 1 


( 10 ) 


Suppose now that from (6) we have 


(k + \)[d(x^ t , x„) + d{x n , x„ +l )]< (2 k + a)d{x n _ v x n ). 


Then we get 


d (x n . x n+l )< ~ a -- d(x nA , x n ). 
k + l 


(id 


Set 


,t =max M±*)±^ J_, ,._L. ,.J±1 

{ h{k + \) Jfc + 1 ifc + 1 J [ h{k + 1) k + l Jfc + lj 


Then by (4), X < 1 and so from (9), (10) and (1 1) in all cases we have 


d{x„, x n+x )<Xd{x n _ x , x n ), 0<X<\. 


( 12 ) 


From (12), by induction, we get d(x„, x n+l )< X n d(x Q , x,). Since 
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Xn)^4 X i’ 

i~m 


* 1+1 


0, 


X,)<A" 


4x 0 ,Xi) 

1-1 


and 1 < 1 , we conclude that {*„} is the Cauchy sequence. Since X is complete, there 
exists some u in X such that 


u = lim x „ . 

oo 


From (5) we conclude that lim,,^ Tx„ = u, too. 
Setting x = u, y = x n in (3) we obtain 


min {hd(Tu, Tx n ), \d(u, Tu)+d(x n , Tx n )]} < (2k + a) 


maxi^(w, x„), —d(u, Tx n ) 

l k 


Taking the limit when n tends to infinity we get 

m\n{hd(j'u, u), d{u, Tu )} < 0 . (13) 

Hence Tu = u and the proof of the theorem is complete. 

Now we shall prove a common fixed point theorem. 

Theorem 2 : Let C be as in Theorem 1 and S', T: C ->C a pair of selfmappings of 
C satisfying the following condition: 


min 


{hd(Sx, Ty), \d{x, Sx)+d(y, 7>)]}< (lk + a)maxyi(x, >>), j-m 



for all x,y in C, where k > 0 is arbitrary and a. and h satisfy (4) and 
m(x,y) = mm{d(x, Ty), d(y, &)}. 
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Then S and Thave a common fixed point. 

Proof : Similarly as in the proof of Theorem 1 we can choose a sequence {jc„} in 
C such that 

d{x 2n . x 2n+x ) = j—^d{x ln , Sx 2n ); d{x 2n+! , Sx 2n ) = y—d(x 2n , Sx 2n ) 

and 

d( x 2n+ 1> X 2 n+l) = 7 , d{x 2n+ \, Tx 2n+] )/ d(x 2n+2 , Tx 2n+\) — ~, 7^(-*2fl+ 1> ^ X 2n+l) 

k + l £ + 1 

for n = 0, 1, 2..., where x Q in C is arbitrary. Replacing x by x 2n+2 and y by x 2 „+\ in the 
condition of Theorem 2, the reader may proceed on the lines of Theorem 1, to derive 
d (x 2 n+i, x 2 n+i) ^ Xd{x ln +\, x 2 n+i), where X is the same as determined in the proof of 
Theorem 1 . Thus, as in the proof of the above Theorem, it can be shown that {x„} is 
the Cauchy sequence and u = lim^^ x n is such that Su = u and Tu = u. 

Remark 1: It is easy to see that a contraction mapping which satisfies (1) also 
satisfies (3) with any a, k and h with a+2k> Xh. In this case assumption of metrically 
convexity of a space is supperfluous. The following example shows that also there are 
mappings which are expansive on some subsets and which satisfy (3). 

Example : 

Let X be the set of reals with usual metric and let T : X —> X be defined as 
follows: 


Tx = -2x, if x > 0 

Tx = -—x, . if x > 0. 

4 

We show that T satisfies (3) on whole X with, for example, a= 1 , k = 2 and h = 6 . 
Let x, y in Xbe such that x<0 and y < 0. Then 
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hd{Tx, Ty) = | d(x, y) < 5 d{x, y) = {2k + a)d{x, y) . 

Let now x < 0 and y > 0. Then 

d{x, Tx)+ d{y, Ty) = ||jc| + 3y < 5^| + y)= {2k + a)d{x, y ) . 

Similarly inequality holds for x > 0 and y < 0. 

Let now jc> 0 and y > 0. Then 

d(x, Tx)+ d(y, Ty) = 3(x + y) . 

Suppose that x<y. Then 

3(x + y) < 5(y - x) = {2k + a)d{x, y), if x<y/4, 

3(x + y)<^-(2x + y)= min{a'(x, Ty), d(y, Tx)\ if y/4<x<y . 

2 k 

Similarly, for x > y we have two inequalities which correspond to the case x>4 y 
and y <x <4 y, respectively. Therefore, T satisfies (3). Since Picard iterations {Tx} 
diverge for any x * 0, we conclude that T does not satisfy (2). T also does not satisfy 
any known contractive condition which implies convergence of Picard iterations. 

Note that any expansive mapping T, defined on X = (-oo+oo) by Tx = -cx with 1 < 
c < 3 , satisfies (3) with a-c- 1, where k > 0 is arbitrary and h satisfies (4). 
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Abstract 

The diffraction of torsional waves by a penny-shaped crack situated in an infinitely long 
isotropic elastic cylinder is considered. Integral transform method is used for finding the 
solution of the equation of motion from the linear theory for homogeneous, isotropic elastic 
material. The problem is reduced to that of solving an integral equation. It is found that only 
torsional waves with frequencies less than a parameter-value, depending in the radius of the 
cylinder, can propagate. The associated integral equation is solved numerically for a range of 
wave frequency and the radius of the cylinder. The solutions are used for calculating the 
dynamic stress intensity factor, displacement on the surface of crack and crack energy. The 
results are presented graphically. 

(Keywords : torsional waves/penny-shaped crack/integral transform/dual integral equation/ 
stress intensity factor/crack energy.) 


Introduction 

Recently great interest has been shown in the problems of the interaction of elastic 
waves with cracks present in elastic solids. The study of these problems is motivated 
by their applications to seismology and exploration geophysics. Elastic waves 
propagating in a solid are modified by the presence of cracks / inclusions and a wave 
diffraction pattern can yield important information regarding the nature of the flow in 
the interior of the solids. 

Due to this reason the problem of finding stress and strain fields in elastic solids 
containing cracks due to oscillatory loading have been extensively studied in recent 
years. Srivastava et al have studied the interaction of shear waves with Griffith cracks 
situated in an infinitely long elastic strip 1,2 and the interaction of longitudinal waves 
with penny-shaped cracks situated in an infinitely long elastic cylinder 3 . An 
interesting conclusion obtained in 1,2 was that only shear waves with frequencies less 
than a parameter-value depending on the width of the wave-guide can propagate in the 
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strip. In this paper we have studied the problem of the interaction of plane harmonic 
incident torsional waves with panny-shaped crack situated in an infinitely long 
homogenous, isotropic elastic cylinder. 

With the help of Fourier sine and Hankel transform of order one a solution of the 
displacement equation appropriate to our problem has been obtained. The boundary 
value problem has been reduced to that of solving dual integral equations, which have 
finally been reduced to a Fredholm integral equation of the second kind. The Kernels 
of the integral equation have been written in the form suitable for numerical 
calculations. As in the case of interaction of shear waves with Griffith crack situated 
in an infinite strip, in this case also, the convergence of kernels of the integrals require 
to impose certain restriction on the frequency of the torsional waves, We have 
concluded that only torsional waves with frequencies less than a parameter value 
depending on the radius of the cylinder can propogate in the cylinder. The solution 
which we have derived are valid only for such frequencies. Numerical calculations 
have been made for a range of wave frequencies and the radius of the cylinder. 
Variation of important physical quantities like stress intensity factor, displacement on 
crack surface and crack energy has been illustrated graphically. 


Formulation of the Problem 

Consider an infinitely long, homogeneous, isotropic, elastic cylinder of radius h, 
containing a penny shaped crack of unit radius whose plane faces are stress free. 
Introduce cylindrical polar coordinates (p,0, z) with the center and axis of the crack as 
origin and z-axis; the crack is then given by z=0,0 < p < 1, and the infinitely long 
Cylinder is given by 0 < p <. h, -oo < z < oo From Love 4 in torsion -situation the only 
non-zero component of displacement is the 0-component which is denoted by U 6 (p,z) 
and the only non-zero stress components are Gq z (p, z) and c p e (p, z) which are given 
in terms of £/ e by 




dU B 

L 

dz 


°pe = PI 


dp p J 


( 1 ) 


where p is the shear modulus of the material. The equation of motion satisfied by U e 
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a 2 E/ e | 1 dU „ 
3p 2 p dp 


U* , d 2 U 6 

p 2 dz 2 


+ k 2 U e =0 


( 2 ) 


where k = eo/c. The speed of shear wave, 



being the density of the 


material. In what follows the time dependence of all field quantities assumed to be of 
the form e m \ the exponential is subsequently omitted in this paper. The plane 
harmonic incident torsional wave in given by the displacement vector (Mai 5 ) 


u° p = 0 = u° 

Ut=B 0 pe' b 


( 3 ) 

( 4 ) 


where B 0 is constant . 

Hence the boundary condition on z =0 are 


U e (p,0) = 0, p > (5) 

dU° 

<*e, (P,0) = -Ce, (p,0) = = ~B 0 p\xik = -q 0 p, 0 < p < 1 (6) 

dz 

If the boundary of the cylinder is stress free, then 


Cf0 2 (h, z) (7) 

The problem of determining stress distribution reduces to that of obtaining 
solution of the displacement equation (2) subject to the above boundary conditions. 

The solution of (2) for an elastic half-space obtained with the help of Hankel 
transform of order one, such that displacement at z = co vanishes is given by 

U e (p, z) = (sgn z) f £/L(£,)e- m J , (£,p )<&, ( 8 ) 


where p 2 £ 2 - A 2 . 
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Again the solution of (2) for an infinite elastic cylinder can be obtained with the 
help of Fourier sine transform. Let 

£/e(p,Q = |* U % (p,z)smt?dz (9) 

By multiplying (2) by sin Cg and integrating with respect to z from 0 to co , we get 

p — ^- + p— (gp+l)E/e=0 (10) 

op ~ op 

where g 2 <£ 2 - A^.The solution of this equation is 

Ue=B(QI i (gp) + C(t;)K l (pg) (11) 

where /j(p) and K\ (p) are modified Bessel functions of the first and third kind. Since 
K\ (p) -»co as p = 0. Hence for a solution for cylinder containing origin C = 0 and the 


required solution is 

U 0 (p,Q = B (C)/,(gp) (12) 

By inversion theorem for Fourier sine transform we have 

t/«(p,2)=|°B(0/|(*p)sm^. (13) 

Thus the complete solution of (2) suitable for the problem stated above is obtained by 
superposing (8) and (13). It is given by 

17 e (p ) z) = (sgnz)f^) e - pi V 1 (^p)^ + f5(Q/ l (gp)sin« (14) 

From stress and strain relation (1) we have 

a 9j (p,z) = -p(sgn z)l^A(^J^P)di + pf ^(Q/,(gp)cos« (15) 
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where 


sgn z = -1, z < 0, 

P = ^ 2 ~k 2 , 

\>,k 


=l,z> 0, 

= -iyJk 2 , 

%>h 

=-/ jk 2 -^ 2 , £>*• 


Derivation of the Integral Equation 
The boundary conditions (5) and (6) are satisfied provided 

f f ^(C)A(gP) = Po p. o < P < l 


p)^ 


l < p < l 


where yf (c,) = c/1© and /? 0 = 7o/p- 

From (7) after using inverse cosine transform, we have 


C5(Q/,(g/0 = l\^A(k) e ^J\mdt\co^zdz 

cos Qzdzflt, 

P 2 


(3 2 +C 




The equation (16) and (17) can be written as 

=P(p), 0<p< 1 

f^(^)J 1 (^p)^=0, P>1 


(16) 

(17) 


(18) 


(19) 

( 20 ) 
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b 2 -k 2 _? 

where H{%) = ^L_ - , H(t) -> as £ oo and 

p( P ) = PoP + f c5(o /,(«p)rfc - f . 


Let the solution of the above equation be 
A'£)= y^q(t)J } fit)dt 

=#f ?(1> X + ^“ si " 5 '4 

With this choice of A'(g) , (20) is satisfied and (19) gives 


(21a) 

(21b) 



1 00 

since ^ (%p) sin %tdt; = 0, fort>p. 


This is an Abel integral equation. Its solution is 




( 22 ) 


By integrating the above equation, we get 
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If we substitute the value of P(p) from (19) in the above equation we get 

V p VP (24) 

-1-7= =fW(WKp)^P ] 

V' 2 -p 2 

By substituting the value of^l'(^) from (21a) in (18) after some manipulation we 


VKOligh) = l^q(t) f dt 

= j£VFg(0 f - C 2 dt 


since C 2 +p 2 =C 2 + ^ 2 -^ = ^ 2 + g 2 . 

By using the results given in Gradshten and Ryzhik 6 [ pp 692, 694 ] we get 

(25) 

By changing the order of integration, (24) can be written as 
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For deriving the above equation we have used the results given in 6 . By 
substituting the values of A' (Q and C ) B{ Q from (21) and (25), then by changing the 
order of integration we get the following integral equation for determining the 
unknown functionx(0- 

XW+ l[L,Xt,m) + L 2 (t,mMt)dt = t 2 (27) 


2 2 

where q(t) = -J-p.x(0, 

3 \ n 

LXt,m) = (28a) 

and Z 2 (t, m ) = ^[^I K (gO/ r Xgm)f^c/(;. (28b) 


Both the functions L\ ( t,m ) and L 2 ( t,m ) are symmetric functions of t and m . 

The integrand in L\ ( t,m ) has no poles. It has only branch point at the point £ = k . 
Following the procedure given in 5 , the infinite integral in L\ ( t,m ) can be converted in 
to the following integral with finite limits. 





= m>t (29) 


Since 




and 


=-/^ 2 -c 2 =-ig\ S <k , g'=ik 2 -tx ' 
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we have 


L 2 (t,m) = 


tm i 
2 

2 


~*g 


~ Ixi-ig'Wxi-ig'm) 




IiHg’h) 


dU 


f — Iy(gt)Iy(gm)^Q-d$ ] 
A g A A Ugh) s J 


tm r fix 


2 1 2 




c 2 


, rt ff, ( v*) 






Mg'h) 


^ + 


tyS&Vy^gm) 


K,(gh) 

Ugh) 


dt; } 


By substituting = A 2 (1-y 2 ) in the first integral and = A 2 (1+y 2 ) in the second 
integral we get 


L 2 (f "0 = J t T I (1 ~ ^ J X ~ -tttt ■■■ - dy 


H?\kyh) . 

UUh) 


+ [(\ + y 2 )' Al 3A (kyt)L A (kym)jj^dy ] 


(30) 


It may be noted that the first integral in (30) is convergent only when the 
dimensionless frequency k is such that kh < j\ \, where y']],=3.83 17060, is the first 
finite zero of J\(x) when the zeros are arranged in ascending order. Hence it can be 
concluded that only the torsional waves with k < yij/h can propagate in an elastic 
cylinder of radius h. This fact is similar to the known result for the propagation of the 
shear waves in an infinitely long elastic strip containing a Griffith crack 1 . 

Quantities of Physical Interest 

By using (15), (21b) and (32) it can be shown that the stress component a 6; .(p,0) 
in the neighborhood of crack is given by 


a ez (p,o) = -n jvgv, gp w + o(i) = + o(i), P > i 
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Defining stress intensity factor by 


N = p Hm Q (P»°)| = 


3n 


(31) 


Another quantity of physical interest is the displacement on the crack surface. 
Using (14), (21a) and (21b) we get 


U 6 ( p,0) = 


4< ?o f PX(0 M _ 4 ?0 ( 
3ttfi -p 2 37tp 

-Af(p), 0<p<l 




4 9o 


(32) 


37tp 


Crack energy can be calculated by using (35). Its value is 

W = 2 nq 0 |pf/ e (p,0)rfp = ^ |V X (0<* = ~ (33) 

Numerical Calculations 

The integral equation (27) has been solved numerically. Using the method of Fox 
and Goodwin 7 (27) has been converted into a system of linear algebraic equation. The 
second integral of Li(t,m) in (30) has been evaluated by using five point Gauss - 
Laguerre quadrature formula while the first integral of L 2 ( t,m ) and the kernel L\ ( t,m ) 
in (29) have been evaluated by using Simpson’s quadrature formula. A complex 
computer programme has been used to evaluate X (0- These values have been used to 
calculate dynamic stress intensity factor; the displacement component U 6 and the 
crack energy. Calculations have been done for a range of dimensionless frequency k . 

Since for k = 0, g = C, the kernel L\ ( t,m ) vanishes while L 2 ( t,m ) from (28b) 
reduces to the following form. 


L 2 (t,m) = 




(34) 
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Two different values of the radius of the cylinder viz. 1.5 and 2.5 have been 
chosen. Keeping the limiting condition k < j-^ylh in view, the values for k 
corresponding to h = 1 .5 are from 0 to 2.5 in steps of 0. 1 . Similarly in the case h - 2.5, 
k has been taken from 0 to 1.5. The graphs shown in Fig: 1 represents the amplitude of 
1x0)1 plotted against the dimensionless frequency k. It is evident from these 
graphs that for these values of the radius h- 1.5 and 2.5, | x(l) I reaches minimum 
value for k = 0.1, then shows an oscillatory tendency. 

In Fig. 2(a), displacement | M( p) | has been plotted, for h= 1.5 against frequencies 
k= 0.4, 1.0, 1.6, 2.0 and 2.4. In Fig. 2(b) displacement curve for h = 2.5 and k= 0.4 
and 1 .4 is plotted. These figures show that the multiple reflections inside the crack 
introduces oscillations in these curves. 

Fig. 3 illustrate the variation of | W(k) | the crack energy with the frequency k. 
For h = 1 .5, I W(k) | first decreases, has minimum value for k = 0.6 then increases and 
reaches maximum value for k - 1 .4, then again decreases with increase in the values 
of the frequency k. While for h = 2.5, | W(k) | first decreases to minimum value for k 
- 0.4 then increases with k. 
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Abstract 

(i) First we consider the structure of the electron. If an electron is a charged body - a sphere, the 
charge on its various parts should repel each other, and the sphere must expand till it blows off. 
To retain its shape, there must be some mechanism. It is the spin. Electrodynamics teaches us 
that a circulating, charge produces a magnetic field along its axis. Now a play between the 
spinning charge and this induced magnetic field may be imagined to provide a dynamic 
equilibrium between the opposing forces and hence the stability of the shape and size of 
electron. Therefore, intrinsically an electron is not only a charge but it is also a tiny magnetic 
dipole as well. For the electron to be dynamically stable it should posses spin angular 
momentum and contain certain amount *of electromagnetic energy and that might be appearing 
as inertial mass. 

(ii) In our view and according to Thomson 1 as well, an electron possesses mass solely by virtue 
of its electric charge and has no mass in the material sense it can be shown that the mass 
possessed by the electron in the frame of reference in which it is at rest gets increased in a 
natural way in the frame of reference in which the electron is moving with an uniform velocity, 
This relation is exactly the same as we get by the relativistic formula for the variation of mass. 

(iii) If we generalize the view that the intrinsic mass is always by virtue of electric charge, as 
every material particle in the universe is composed of nuclear and electronic charges, it should 
vary with the velocity of the particle. Therefore it suggests that the mutual gravitational 
attraction of the masses must also be a kind of electromagnetic interaction since masses too 
have some kind of electrical origin, This view provides us a new way to picturise the planetary 
motion as well. Possibly it provides the answer to the problem that why is the elliptical path is 
acquired in general by a planet. 

(iv) Now for a gas. a molecule is an entity which is blanketed by a dynamic electron in the 
outermost orbital of the molecule as a whole, The electronic charges, their orbital and spin 
angular momenta and the nuclei satisfies the valencies and keep up the molecule electrically 
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neutral in the dynamic equilibrium state, Except the outermost dynamic electron which 
envelopes the nuclei with the other satellite electrons within the molecule, offers itself as an 
enclosure to a!! confined nuclei and electrons and completely shields their electromagnetic 
interactions out side the enclosure and represents the only dynamic electron capable of 
influencing the similar electron of neighbour molecule by inertial interaction. Thus in a gas, 
every entity may be regarded as a mass Point of the rest mass, mass equivalent to that of the 
electromagnetic or electrical mass of an electron providing some sort of very weak, attractive 
force of electromagnetic or of electrical origin acting on its similar neighbor molecule / entity, 
which may not be other than the well known gravitational force. 

(v) Along with this we may also consider one gm mole of a gas as a single whole of a bulk of 
material at a given temperature, which is under thermal equilibrium with the surrounding 
radiant energy in which the photons of, various frequencies following the temperature 
dependent statistical laws of distribution are under motion in all possible directions. 

(vi) In this case there is very large empty space between the entities which may be regarded as 
free space provided for the motion of the photons and when they have occasional encounters 
with the outermost electrons, photons get absorbed by them followed by their release and so on. 
These consideration yield the formula for Universal Gravitational Constant, G in terms of 
standard values of velocity of light, rest mass of an electron, Stefan's constant; Avagadro's 
number, the volume and temperature of one gm mole of a gas at S.T.P. which is given bellow : 

G = 21.493226 v.T*V 0 4n / N* n .m^c 

The values of G when calculated comes out to be 6.6627 X 10‘ 8 c.g.s. which is strikingly v ery 
close to the standard value of gravitational constant G and just falls lower by 0.157 % only 
supporting its origin in the assumptions made earlier. 

(Keywords : relativity/planetary motion/molecular gas/gravitation/electron). 


1. Introduction 

The Possible Structure of an Electron : 

The set of equations of Lorenz transformation was used by Larmor to explain the 
null result of the Michelson - Morley experiment, in his Aether and Matter, pp 174 - 
176 (Cambridge University Press, New York 1900). Later the out come of the 
Michelson - Morley experiment as a new principle as that the speed of light 
(electromagnetic waves) is independent of the motion of. the light source or receiver 
formed the basis for the search of the equations of transformations for a frame of 
reference moving uniformly relative to other frame of reference. This lead scientists to 
consider it m two ways, (i) Lorenz considered a light wave spreading out from a point 
source at the origin in a reference frame. The space being attributed with the 
properties i.e.. it is isotropic and homogeneous in nature, the wavefront (surface of 
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qual phase) will be a sphere. But according to the new principle the wavefront must 
also be a sphere when viewed in an another reference, frame which is in uniform 
motion with respect to the source, (ii) While as the bodies are perceived in Nature by 
light signals Einstein tried to resolve the problem, adhering to the principle of 
equivalence and the constancy of the velocity of light resulting to the same Lorenz 
equations of transformations for a frame of reference moving uniformly relative to the 
other frame of reference. These equations do not concern with the intrinsic nature of 
the objects. Even an electron which has charge and the spin both-how does vary in 
mass while it is under relative motion? What happens to its structure? These are the 
questions which remain unsolved. If an electron is a charged body - a sphere, the 
charge on its various parts should repel each other and the sphere must expand till it 
blows off. To retain its shape, there must be some mechanism. It is the spin. 
Electrodynamics teaches us that a circulating charge produces a magnetic field along 
its, axis. Now 'a play between the spinning charge and this induced magnetic field 
may be imagined as to provide a dynamic equilibrium between the opposing forces 
and hence the stability of the shape and size. 



Fig. 1 
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2. The Spin Angular Momentum, the Mass of the Electron and 
its Variation with Velocity 

Consider an electron of mass m 0 and as an electronic charge e smeared over on a 
spherical surface of radius r e and say spinning with the average peripheral velocity v, 
= 27t r e / T s about X - axis as shown in the Fig. 1. The de Broglie wavelength X 
associated with such an electron will be given by X = h/moV s , where if X/2 be equal to 
2nr e , then obviously spin angular momentum will be given by 

mov s r e =h/4n = s.(h/2n) 

where s = 1/2, a result used to detail the atomic and molecular physics. This spin 
motion of time period T s will generate a magnetic field 

H= fx 0 0 / T s ) / 2 r e = ( \x 0 ev s /4nr 2 ) 

along the A- axis which in turn (due to electromagnetic interactions of magnetic field 
with the charge will keep the electronic spin motion sustained. In this case the Lorenz 
force will act as the centripetal force equal to Hev s which in turn must be equal to 
centrifugal force moV s 2 / r e . Equating these two we get 

Hev s = m 0 v ] jr e or (p 0 ev s /4nr 2 ) ev s = m Q Vi 2 /r e . 

substituting po = 1 / s 0 c 2 , we get that the electrostatic energy, 

E = e 2 /4n s 0 r e = m^c 2 

Obviously for this model the electronic mass appears as the equivalent mass to its 
net electrostatic energy, concluding that the electron the ultimate particle of 
electricity, has no mass in the material sense. Its mass is purely electrical. 

Now let an electronic charge be smeared on a spherical surface of radius r e and 
say spinning about X-axis as shown in the Fig. 1. Let O', the origin of S frame of 
reference, be the centre as well of the charged spherical shell which is moving with 
uniform velocity v along X-axis with respect to another S frame of reference having 
the respective axes parallel to each other and being its origin O is just coinciding with 
O' at the instant the figure has been drawn. This three dimensional spherical surface 
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may be generated by revolving a circle of radius r e in the X'-Y' plane about the X-axis 
with its centre at O'. The equation of this plane of intersection of the spherical surface 
with the X- 7'plane will be given by 

x 2 +y 2 =r? (1) 

Let in the X-7'plane, P'{x’. y', 0) be a point on the charged spherical shall having 
its foot of the perpendicular on the x-axis at C. The circulating charge at P 'and the 
magnetic field acting along the x-axis must be in conformity with each other, that is in 
the dynamic equilibrium the disturbances caused by them will travel with the velocity 
of light c and will be received in pace. Let S t be the time elapse in which a signal 
reaches from C' to P'then 

8 1 = { C'Pyc 

or y'= C'P'= c. dt = r e sinG' (2) 

where x'=r e cos 9' 

As this S' frame of reference is in relative uniform motion with S' frame of 
reference the position of C'will appear at S further away from c by a distance 

CC' = (v St) = {C'P') vie = r e sin 6'. v/c (3) 

Now due to invariance of the speed of the light in uniformly moving frame of 
references with respect to one another, a result which is well confirmed by Michelson 
- Morley experiment, the point P 'will appear at P so that . 

CP = CP' = y' = r e sin0'. (4) 

Therefore the ordinate y of the point P will appear into the S frame of reference of 
length CP and will be given by 

y = {{CP) 2 - {C'C) 2 ) V2 = {{r e sinQ') 2 - {r e sin 0'.v/c) 2 } 1/2 
= r e sin0' {1 -(v/c) 2 } 172 
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= /{l -(v/c) 2 } 1 ' 2 (5) 

while the abscissa x of P will-remain the same as that of F i.e. 

x = OC' = x /=z O'C r e cos©' (6) 

The equation of the plane of intersection, of the Surface with the X-Y plane as 
observed in the S frame of reference will be had by substituting the values of. x' and/ 
in terms of x and y into the eqn. (1) then we have 

x 2 +y 2 {l-(v/c) 2 }’ /2 = r 2 
or x 2 /[r? ] + y 2 / [r 2 {l - (v/c) 2 } } = 1 
or x 2 la 2 + y 2 /b 2 = 1 

where a = r e and b-r e {l-(v/c) 2 } 1/2 . 

When such an ellipse is rotated about x-axis it will generate an ellipsoid having its 
minor axis perpendicular to the direction of the motion. This gives us the structural 
change of a spherical shell of electronic charge which is at rest to an ellipsoidal shell 
of the same charge moving with an uniform velocity v. Thus the electronic charge 
smeared over the curved surface of the area 4na~ of sphere of S’ frame of reference 
will now appear as if it has smeared all over the curved surface of the area 4n ab of 
the ellipsoid in S frame of reference. Naturally the concentration of the charge will 
increase in the parts of the same quantum of charge, i.e. in the parts of the electron. 
This compression of the charge with in the body of the electron should increase the 
electrical energy stored in the surface possessed by the electron itself by the ratio of 
the curved surface area of the sphere to that of the ellipsoid. As the curved surfaces of 
the sphere and the ellipsoid are 4 na 2 and 4 nab respectively. Therefore the ratio of the 
area of the sphere to that of the ellipsoid will be alb. As the electric energy of the 
sphere "of radius a is e 2 /47tSotf, the energy of the ellipsoid may be estimated as given 
below : 

The energy of the ellipsoid = (e 2 /4m 0 a) alb = e 2 /4xe 0 i 

= The electrical energy of a sphere of radius b. 
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Thus the electrical energy of a charged ellipsoid numerically is equal to e 2 /4ns 0 b 
which is just the same as that of a charged sphere of radius b. Let us assume that the 
charge spherical shell of mass m 0 radius a being in the S frame of reference and the 
charged elliptical shell of mass m v in the S frame of reference, then their electrical 
energies will be 

E a = e 2 /4Tcs 0 a and Eb = e 2 /4mo b respectively 
The ratio E b : E a will be given by 
E b !E a = alb = 1/d-vV) 1 ' 2 , 


where a = r e and b = r e { \-v/c) 2 } 112 . 

Further as an electron possesses mass solely by virtue of its electric charge and 
has no mass in the material sense therefore the ratio E b : E a must be equal to the ratio 
m v : m 0 , thus 


m v lm 0 = E b /E a = a/b = l/( 1— v 2 /c 2 ) ,/2 


or m v = « 0 /(l-v 2 /c 2 ) 1,J . 

This electron structure based formula , tells us that the mass possessed by the 
electron in the frame of reference in which it is at, rest gets increased in a natural way 
in the frame of reference in which the electron is moving with an uniform velocity but 
so for the contraction in the radius of the charge sphere is concerned, it is 
perpendicular to the direction of motion. The theory of relativity derives the same 
relation for the variation of mass. In the relativistic formula the contraction takes place 
only in the direction of motion as in this case, a sphere is carried by an observer in 
motion it will appear as an ellipsoid for an observer at rest with its minor axis along 
the direction of the motion, where the structure of the particle is in no way considered. 

3. The Planetary Elliptical Motion 

As we have observed that the mass is purely electrical in nature and when this 
charge is in motion its mass increases. Every particle in the universe is composed of 
nuclear and of electronic charges having a mass and therefore it should vary with the 
velocity of the particle. This variation in mass must be numerical^ equal to that would 
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have obtained due to the application of the spatial theory of relativity. Further it 
suggests that the mutual gravitational attraction of the masses must have some kind of 
electromagnetic interaction since the mass itself have some kind of electrical origin. 
This view provides us a new way to picture the planetary motion. 

A planet of mass m p , in general moves on an ellipse around the Sun of mass M s 
being at the focus Now whether the planet moves on a circular or an elliptical orbit the 
inverse square law of force will hold. The circular motion is more symmetric than the 
eccentric elliptical motion. It suggests that there must be one more condition which 
enables a circular motion to adopt the elliptical motion. In general a celestial body 
with initial velocity V ps , with respect to Sun might have come from very remote place 
of the galaxy and entered into the gravitational field of the Sun which later on has 
appeared as a planet of mass m p . 

Let us consider first that a planet of mass m p is moving on a circular orbit of 
radius around the Sun of mass M s . For the dynamic equilibrium the centripetal force 
( GMsiripla 2 ) must be equal to the centrifugal force (m p V 2 /a) and the equation of the 
path will be 

x 2 + y 2 = a 2 . 

If this has initial velocity Vps with respect to Sun, the Sun will also has a relative 
velocity V sp with respect to the planet equal and opposite to the velocity V ps of the 
planet with respect to Sun. As we have considered that the gravitational interaction 
also being of some form of electromagnetic interaction, so every action acted by one 
on the other will also travel with the velocity of light c. The relative velocity V sp of the 
Sun with respect to the planet will modify the equation of the path of the planet 
similar to the that which we have observed in the case of the electron. Considering the 
Sun at the origin let x\ / be the coordinates of the planet. Due to the relative motion 
the ordinate y will, get contracted to y (1 -V 2 lc 2 ) m and will become y' while x will 
become x' that is 


as 


x = x' and y = y'l (1 -V 2 /c 2 ) U2 

Substituting these values into the above-equation we get the equation of the path 
x a l[a 2 } +y' 2 l[a 2 {\-(V sp lc) 2 }] = 1, or x l2 /a 2 + y' 2 /b 2 = 1, 
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which is an equation of ellipse with the major axis as a and the minor axis as b = a 
{\-{Vsplcf} m The eccentricity e will be given by s = V sp fc. Thus we find that the, 
orbits of the planets in general are elliptic where sun occupies a focus which is 
directed by the respective relative velocity V sp . 

As both the Sun and the planet are moving with velocity | V sp | relative to each 
other their respective masses will become MJ(\-V 2 p lc 2 ) m and mpl(\-V* p lc 2 ) m let at r 
distance apart, they settle in the mutual fields so that the planet continues to revolve 
round the Sun. Then the centripetal force must be equal to 

G{M s /(\-V 2 /c 2 ) U2 } {m p /(\-V*/c 2 ) ]l2 }/r 2 = [GM s m p /r 2 ] (1 + V*/c 2 ) 

=[GM s m p /r 2 ] (1+e 2 ). 

When a body of mass m p brought with infinitely slow motion from infinity to a 
distance r far from the Sun of mass M s , there is a decrease in potential energy by 
[GMsirip/r], The initial mass equivalent to the potential energy will be A M PE = 
[{ GM s mp/r}/c 2 ]. 

The gravitational interactions have settled the body as a planet of mass m p moving 
with an orbital velocity V p , while in the, initial conditions the body was moving with 
the velocity V ps hence the inertial mass equivalent to the potential energy A Mp# = 
[{GM s mplr}/c 2 ] will also continue orbiting with the same velocity V ps with the planet. 
Allowing these considerations one can calculate the net centrifugal force as 

[m p ] V 2 /r + [A M p ,e] V p Jr = m p V;fr + [{GM s m p /r}/c 2 ] V 2 Jr 


= m p V p /r + [ GM s m p lr 2 ] ( V*/c 2 ) 
= m p V*/r + [GM s m p /r 2 ]c 2 


For equilibrium the centripetal force must be. equal to the centrifugal force, 
therefore 


[GM s m p /r 2 ] (1+e 2 ) = m p V 2 /r + [ GM s m p lr 2 ] e 2 
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or GMstrip/r 2 = m p V* tr 

or V; = GMJr. 

This equation shows that once the body is trapped as a planate: its motion can be 
expressed in terms of the rest mass of the sun and the distance. The only clue 
regarding its initial relative velocity V ps , which is just before appearing the scene, is 
given by tile eccentricity s of the elliptical orbits. The elliptical orbit has more energy 
over and above than that of a circular orbits of radius a (the semi major axis ) by a 

value AMp'E = F^/2 where A is the inertial mass equivalent to the potential 

energy. Similarly the rosette shaped path of mercury about the sun that is the 
processional motion of the major axis of the orbit may be accounted by the variation 
of the inertial mass equivalent to its kinetic energy with its orbital velocity since the 
change in mass is greatest at the perihelion, it results in a precession of the orbit of the 
planet. The perihelion of the ellipse rotates as the time passes on i.e, the major axis of 
the ellipse turns about the sun in the plane of the ellipse and the resultant path of the 
planet is not a simple ellipse but a complicated curve known as rosette. 

4. The Molecular Gas 

Let us consider one gm mole of a gas as a single whole of a bulk of material at a 
given temperature which is under thermal equilibrium with the radiant energy all 
around specified to this temperature. Further an entity is defined by the volume 
v 0 =F 0 /A^ i where the volume V 0 at the standard temperature T 0 and pressure p 0 contains 
N a the Avagadro's number of molecules. From the knowledge of quantum theory a 
molecule is an entity which has blanketed by a whirling electron in the outermost 
orbital of the molecule as a whole. The electronic charges, their orbital and spin- 
angular moment in a molecule with the nuclei keep up it electrically neutral and 
satisfy the valances in the dynamic equilibrium states. The outermost dynamic 
electron which envelopes the nuclei with their other satellite electrons of the molecule 
offers itself as an enclosure to all confined nuclei and electrons as well as it 
completely shields their electromagnetic interactions outside the enclosure except the 
only outermost election is capable to influence the neighbor. Thus every molecule as 
an entity may be regarded as a mass point of the rest mass m e equivalent to that of the 
electromagnetic or electrical mass of an electron providing some sort of very weak 
attractive force of electromagnetic or electrical origin on its similar neighbor 



GRAVITATIONAL FORCE AS THE INERTIAL INTERACTION 


383 


molecule/entity. Let thjs arisen force be the well known gravitational force To effect 
this for a gas at any temperature T the following two assumptions are made viz. 

(i) that ( T/To ) times an entity is the number of entities for. electromagnetic mass 
interaction. 

(ii) that the number of entities present in a gas at any temperature T are defined 
by the ratio V t /Vq, i.e. Vt/V 0 = N A (T/T A ) entities in all. 

At the same time these entities may be considered as surrounded by the space, just 
like that of a black body, which is full of radiant energy in which the photons of 
various frequencies following the temperature dependent statistical laws of 
distribution are under motion in all possible directions. In our case there is very large 
empty space between and within the-entities as molecules being of very small, 
dimensions compared to that occupied by an entity, which may be regarded as free 
space provided for the motion of the photons and when they have occasionally 
encounters with the electrons get absorbed by the electrons of the outer most orbital 
and then followed-by. their release and so on. Consider the ideal gas at S.T.P. having 
Avagadro's number N A (6.022169 x 10 23 ) of entities say primary molecules in volume 
V 0 (22410.656 cc) at temperature T 0 (273.16 K). Further we may estimate the space 
volume Vo provided for this single entity. For to evaluate this simply the volume has 
to be divided by the number of the molecules/entities present i.e. V 0 = V 0 /N A = 
3.7213595 x 10~ 20 cc per entity be a constant quantity strictly for an entity at all 
temperatures. Let a 0 be the radius of average spherical confinement for this entity of 
volume equivalent to Vo then = 4n a\ / 3 , therefore ao — 20.71 1156 x 10 8 cm. 

As the size of a molecule is a few Angstrom (10' 8 cm), the space in which an 
entity lies is thousand times larger than that of it. Therefore, one may picture this 
situation classically as eveiy entity with regard the other entity is situated very far 
from one. another. The outermost dynamic electron which envelopes the nuclei with 
their other satellite electrons of a molecule, offers itself as an enclosure to all confined 
nuclei and the other electrons as well as it completely shields their electromagnetic 
interactions outside the enclosure. The only whirling electron is capable to in 
influence the neighbor, every entity may be regarded as a mass point (of the rest mass 
the equivalent to that of the electromagnetic mass of an electron) providing a some 
sort of very weak inertial attractive force on its neighbor molecule/entity arisen due to 
some kind of unknown electromagnetic interaction with similar entities at :S.T.P. Let 
this force be the well known universal gravitational force. 
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5. Tfae-Molecule and Photon Interaction 

The study of blackbody radiation has revealed that the space inside the black body 
is full of radiant energy in which the photons of, various frequencies following the 
temperature dependent statistical laws of distribution are under motion in all possible 
directions. In our case there is very large empty space between and within the entitles 
which may be regarded as free space provided for the motion of the photons and when 
they have occasionally encounters with the electrons get absorbed by the electrons of 
the outer most orbital and then followed by their release and so on : 

6. Equilibrium between the Thermal Radiation and Potential 
Energy of Molecules 

Combining the concepts of [4] and [5] together we may consider one gm mole of 
a gas as a single whole of a bulk of material at a given temperature which is under 
thermal equilibrium with the radiant energy all around specified to this temperature 

To resolve this problem — 

(a) we first calculate the lowering of potential energy of an entity as being under 
the action of a weak inertia! attractive forces viz. universal gravitational 
forces of other entities 

(b) Secondly we estimate the number of entities enveloping the appropriate shell 
which is holding one gm mole bulk. 

(c) The net potential energy of these entities, which may be regarded as the 
ground state of the bulk, may be excited by the absorption of the incident 
radiant photons and be re-radiated out after the relaxation time. This released 
energy per unit cross section of the exposing surface of the bulk per unite 
time must be equal to the amount desired by Stefan's law for thermal 
equilibrium. 

6 (a) The lowering of potential energy of an entity as being under the action of a weak 
attractive gravitational influences of the rest of entities 

First we consider a gas at S.T.P. The volume V 0 at the standard temperature To and 
pressure Pq contains .N A number of entities. The volume Vq = .Vo / N A may be 
considered as a volume of a cell which is necessary to define one entity at 
temperature T 0 . Now let us introduce two assumptions as follows : 
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(i) that T / Tq times an entity is the number of entities for electromagnetic mass 
interaction, 

(ii) that the number of entities present in a gas at any temperature are defined by 
the ratio Vjl To, i.e. Vf! To = N A (T / To) entities in all. 

Let cit be radius of the spherical shell provided for a single molecule at 
temperature T then V T = V 0 T / T 0 = N A V 0 T / N A T 0 = (N^T/T 0 )(4rr/3)a/ a \ . where a 0 = 
(3 JV4rt A^) !/3 . Further let us draw imaginary concentric spheres of "radii a 0 , 3 a 0 , 5 

ao, (2r-l) ao to (2n r \)ao, so that the last sphere of radius (2«r 1) tfo, 

includes the whole volume of 1 gm mole of the gas at temperature T, then N A number 
of molecules which were the number of entities at temperature To will now assume the 
number N A (T/T 0 ) then. (4rc/3) [(2« r l) a a ,] 3 = (N A T/T Q ).(4n/3) a\, or (2n r \) = 
(N A T/To) m . N a being a large quantity. 

2 n T - (A^r/7o)' /3 + 1 and will be approximately equal to ~ (A^ 777V) 1/3 or n T = (1/2) 
(N A T/T 0 ) m =42223509 (77T 0 ) I/3 . 

Thus the radii a 0 , 3 a 0 , 5 a 0 , , (2n r l) a 0 will divide the whole space in 

(n r 1) concentric spherical annular zones z 0 , z u z 2 , etc. in which the rth zone (z r ) will 
be comprised of {2r+l) 3 - (2r-l) 3 } entities. These entities will be equidistant from the 
centre of the central zone Zo and the centre of this zone z r will be at a distance r.2ao 
from the centre. 

Now to visualize a spherical shell filled with one gm mole of the gas, we first 
assume that initially there is a single entity at the centre, then z\ zone of which mid 
point is at 2 a 0 , and continuing to say to rth zone z, filling it with.{(2H-l) 3 -(2r- 
l) 3 } entities which will be at 2rao distance till r becomes wy-for the last zone. 

To calculate the lowering of potential energy of an entity of this system, we 
assume first that up to z r .\ th zone have already been filled with (2r-l) 3 entities. Using 
the assumption (i) that (77 T 0 ) times an entity is a factor for electromagnetic inertial 
interaction and the only the rest mass of the outermost whirling electron of an entity is 
capable to influence the neighbor, the effective inertial mass of this core will become 
[(2M) 3 . 77 To mo}. Further to fill the z r th zone, we have to bring {(2r+l) 3 -(2r-l) 3 } 
number of entities, in which every entity will be at a distance 2rao from the centre. 
The effective inertial mass of this core will become 

[{(2r+l) 3 -(2r-l) 3 }(7’ / T 0 ). w 0 ] . • 

The drop in potential energy A W r will be given by 
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AW r =G[2r-\fT/T 0 . m 0 ] {{2r+\f-{2r-\f}TIT Q mo\l2ra 0 , 

= Gm\{T l To) 2 (2r- l) 3 . (12r 2 + \)/a 0 .r 

= G[ m\ ( T / To) 2 / fl 0 ] [96/ - 144 / + 80/ -24/- +1-1 / r] 

Therefore the net potential energy of the entity at the center may be had by 
summing up the above expression for r = 1 to n T , i.e. we have 


W = ^AW r = G[ ml (T/ T 0 ) 2 /a 0 ] f)[96r 4 - 144r 3 +80r 2 -24r+l-l/r] 

r-\ r= 1 

= G?« 0 2 (r/ro) 2 /ao.] [96/5] « 2 

neglecting cubic and lower powers of r since n T = [1/2] (( N A T/To) m ) being a large 
quantity and a 0 = (3 V 0 /4N A ) 113 Therefore we have 

fT= [ 1/2 5 ] [96/5] [GN 2 A m 2 0 (T/ r 0 ) 1 1/3 ]/(3 F 0 /47t) I/3 

The value of W so obtained for single entity may be generalized to every entity all 
being under the dynamic equilibrium. 

6(b) Estimation of the number of entities per unit area of the enveloping surface of the 
appropriate shell which is holding 1 gm mole of the gas : 

For the estimation of this number we divide the whole volume V T of the spherical 
shell into tiny cubes having their faces normal to X, Y, Z axes. Let 5 be the edge of 
cubical cell of volume V 0 at the temperature To then Vo / N A = S 3 , which confines a 
single molecule. Now as we know that 


V T = VoT/To = N a V 0 T/N a To = (N a T / T 0 ) (V 0 /N A ) 

= (N a T/To) 5 3 = (4jt/3)[(3/4jt)(JV^/8) (T/T 0 )] (85 3 ) 
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where n* = [(3/4n) (N A I 8 ) (77 r 0 )] 1/3 

This gives us a clue and proposes how to construct a spherical shell of radius 
n*.2§ having its volume V T made by tiny cubical cells of edge 28 arranged along the 
axes X, Y, Z. In this case we can calculate the number of the cells constituting the 
peripheral enveloping spherical surface similar to previous section and will be equal 
to 

{(2m* + l) 3 - (2m*- l) 3 } = 24 m * 2 + 2 = 24 n* 2 approximately 

Using the assumption (i) that (T / T 0 ) times an entity is a factor for electromagnetic 
inertial interaction, the effective number of the cells constituting the peripheral 
enveloping spherical surface will become .24m* 2 (T / 7o).These peripheral enveloping 
cells being in thermal equilibrium with the outer radiation, they can casually absorb 
and permit the incident radiation to pass through to the interior to get absorbed by the 
other molecules during their passage through the bulk ,by which the gravitational state 
of these entities may get excited. As the exposed surface is 6( 25) 2 of a tiny cell and 
N a Tl To, according to the assumption (ii), is the number of entities all over the 
volume, the net surface will be equal to 6N A (T / To) (28) 2 The average number of 
entities per unit area of the exposed surface 

= 24 m * 2 (77 T 0 )/6N a (77 T 0 ) (25) 2 = m * 2 / A^5 2 . 

Therefore the potential energy of these entities is lowered by W (n 2 IN A 8 2 ) per 
unit area normal to any of the axes where m* = [(3/4rc) (N A /8 ) (T / 7o)] 1/3 , 8 = (Vo/ 
N A ) m and W=[ 1/2 5 ] [96/5] [GN 2 m 2 0 (T/T 0 ) UB ]/(3V 0 /4n) m . 

6 (c) The relaxation time A t and deduction of the formula for G : 

This net potential energy per unit area normal to any of .the axes due to attractive 
actions of the entities, which restrict the free motion of the molecules, may be over 
come by absorbing the radiation in its course of passage through a cell and re- 
radiating out it after the relaxation, time At to maintain the thermal equilibrium. The 
relaxation time may be considered as the time required to travel a photon though the 
distance 28 r where 28 r = 2(V 0 /N A f 3 (T/T 0 ) m being the edge of an actual cubic cell at 
temperature T therefore At - 28 T /c . 
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Therefore the energy released per unit cross-section per sec. will be given by 

(n‘ 2 /N A 8 2 ) W(c/25 r )=( 1 IT) n 2 cWIN A b 2 8 r =(1/2) n 2 cW/V 0 ( T/T 0 ) u 3 

This released radiation energy amounting equal to potential energy per unit area 
per sec normal to any of the axes of 1 gm mole of the bulk must be equal to the 
amount received by the Stefan's law ( the normal radiation per writ area per sec = 
a I 4 ) for thermal equilibrium. Therefore equating these two we get 

a T 4 = (1/2) n 2 cWIVo (777o) 1/3 

= (1/2) [3/471^/8) (77 T 0 ) 2/3 .c.[1/2 5 ] [96/5] 

[ GN]m 2 0 {T/T)" n ]/(3 F 0 /47t) ,/3 . V 0 (T/T 0 ) ]n 

or a = G [{3/4xN A } M /8V 0 ] [1/2 5 ] [96/5] N 2 A m 2 (l / r o ) 4 ]/(3V 0 /47t) 1/3 , 
or G = a/(3Fo/47t)'7 [{(3/47t)^} 2/3 /8F 0 ]c [1/2 5 ] [96/5] iV> 0 2 (l/r 0 ) 4 ] 
or G = 2 1 .493226 a. T 4 V 4n / N* A n m\ .c 

The value of G may be calculated using the standard data 2 given below : 


c - 2.997925 x 10 10 cm/sec. 

Lf 

mo = 9.1095 x 10~ 28 gm 

M 

N a = 6.022169 x 10 23 



V 0 = 22410.656 cc ]} 

To = 273.16 K K 

ct = 5. 6697 x 1(T 5 ergs per sec. per K 4 


G = 6.6732 x KT 1 


c.g.s 


M~ l L 2 T 2 
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The value of G when calculated comes out to be 6.6627 x 10~ 8 c.g.s. which is 
strikingly very close to the standard-value of gravitational constant G and just fails 
lower by 0.157 % only supporting its origin in the assumptions made earlier. 


Conclusion 

Considering the above said the mathematical calculations yield the formula for 
universal gravitational constant G in terms of standard values of velocity of light, rest 
mass of an electron, Stefan' constant and Avagadro's number of molecules, the volume 
and temperature of one gm mole of a gas at S.T.P. which as given below : 

G = 21.493226 a T 0 4 ,V 0 An / N* A n .m 2 0 .c 

The result is in fair agreement with the standard value of G. 
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Abstract 

The dependence of the values of rotational constants B e of the diatomic oxide and chalcogenide 
molecules of the elements of the groups IIA, IIIA, IV A, IVB and VA on electronegativity has 
been studied. The equation proposed by us earlier for the diatomic halides, in which the bonds 
exhibit predominant ionic character, is found to be valid for the IIA and IVB oxides and 
chalcogenides. Tor the rest of the oxides and chalcogenides, the equation proposed by us for the 
intragroup and homonuclear diatomic molecules is applicable indicating that the 
electronegativity equalization takes place in these molecules giving a covalent character to the 
bonds. 

(Keywords : rotational constants/electronegativity/diatomic oxides/chalcogenides). 

Introduction 

The study of molecular spectra provides very useful information about the 
electronic structure as well as the vibration and rotation of the atoms in molecules 1 . 
Such studies in turn help in obtaining accurate information about geometrical 
arrangement of atoms in molecules, leading to the understanding of many of their 
physical and chemical properties. 

Several properties of molecules depend upon the electronegativity values of their 
constituent atoms 2 " 6 . The equilibrium intemuclear distance r e in a diatomic molecule is 
decided by the force of attraction between the two atoms, which can be considered to 
be of related to the electronegatvities of the constituent atoms. The intemuclear 
distance r e is also related to the rotational constant B e through the expression 7 
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where h is the Planck's constant, c is the speed of light and p. is the reduced mass of 
the diatomic molecule. In our earlier papers we have attempted to systematise the data 
on rotational constant B„ of some diatomic halides 8 ' 9 and homonuclear and intragroup 
heteronuclear diatomic molecules’ 0 with the help of electronegativity. We have also 
been able to estimate the values of B e in those cases in which the spectroscopic data 
on B e were not available in literature. In the present paper, we have made an attempt 
to see if the available data on B e in diatomic oxides and chalcogenides can be 
systematised in a similar fashion. We have considered only those molecules here, on 
which some work has been reported in literature, and which, therefore, seem to exist. 

Rotational Constants and Electronegativity 

Our earlier work shows that the rotational constant B e of molecules exhibiting 
predominant ionic character can be correlated with the electronegativity difference Ax 
of the constituent atoms A and B through the expression 8 ' 9 


log B e = m log 


1 ' 

vp 2 ; 


+c. 


( 2 ) 


We have tested this expression for the diatomic halides of the elements of the 
groups IA, IB, IIA, IIIA, IIIB, IVA and VA, in which the halogen atom is an electron 
acceptor. The values of the slope m and intercept C for the plots of log B e versus log 
(Ay/p 2 ) are found to be different for each of the groups. 

For the intragroup and homonuclear diatomic molecules however, the situation is 
different. In the case of a homonuclear diatomic molecule AA, since the 
electronegativity difference is zero, equation (2) is not applicable. We have shown, 10 
that in the case of such molecules of group IA, VA, VIA and VILA, the following 
equation 


log B e = m log 


v 

vp 2 ; 


+c, 


( 3 ) 


is applicable. With the help of this equation we have estimated 10 the data on B e in 
some homonuclear diatomic molecules, viz. Rb 2 and Te 2 which were not reported in 
literature". 
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The intragroup diatomic, molecules are strongly covalent because the constituent 
atoms have a strong tendency of sharing electrons rather than transferring them from 
the element of lower electronegativity value to the element with higher 
electronegativity value. In the case of a diatomic molecule held together by a single 
covalent bond, if the constituent atoms are different, then the polarity to the covalent 
bond is imparted to an extent and direction indicated by the initial electronegativity 
difference A%. Sanderson 12 has suggested that to form a covalent bond, this must stop, 
and the electronegativities of the two atoms eventually become equal. Our earlier 
work has shown 10 that in intragroup heteronuclear diatomic molecules AB, the 
rotational constant B e is related to the electronegativity values Xa and Xb of the 
constituent atoms through the relation 


log B e = m log 



\ 



J 


+C. 


( 4 ) 


It can be seen that equation (3) is just a special case of equation (4). 


Rotational Constants of Diatomic Oxides and Chalcogenides, 
and Electronegativity 

In the present paper, an effort is being made to see whether equations (2) and (4) 
can be extended to the diatomic molecules in which the acceptor atom B can accept 
two electrons. The diatomic oxide and chalcogenide molecules are expected to exhibit 
such behavior, and hence the present study has been undertaken. As a spin-off of this 
study, we have been able to estimate the values of B e . for MgS, SrS, SrSe, TiS, InO, 
AlSe, GeSe, SiTe and SbO, for which data are not available in literature. 

In Table 1 are given the values of the reduced mass p as cited by Huber and 
Herzberg 13 and of the rotational constant B e of the oxides and chalcogenides of the 
groups IIA, IIIA, IV A, IVB and VA, as cited in the CRC Handbook of Chemistry and 
Physics 11 , 80 th edition (1999-2000). The data on B e on the molecules BiO, BiS, PS and 
HfO which are not given in the CRC Handbook are taken from reference 13. These 
values are reported to be either uncertain or of B 0 , the lowest vibrational states, as 
cited by Huber and Herzberg. The data on the oxides and chalcogenides of other 
groups are too scanty. In this Table are also given the electronegativity values Xa an ^ 
Xb (which are taken from the scale of Mande and coworkers 14 ’ 15 , as was done in our 
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earlier papers 7 ’ 10 ) along with the values of Ax and ^XaXb • have plotted the 
values of B e against the values of log (Ax/p 2 ) for IIA, IIIA, IV A, IVB and VA oxides 
and chalcogenides. We obtain linear plots for IIA and IVB oxides and chalcogenides 
(Figs land 2). However, the scatter in the plots for IIIA, IVA and VA chalcogenides is 
found to be large. Hence these curves are not shown in this paper. 


Table 1-The values of electronegativity, reduced mass \x and rotational constant B e of different 
molecules. 


Molecules 

lA 

1b 

Ax 

JXaXb 

|u(amu) 

B e (cm'') 

IIA oxides and chalcogenides 






BeO 

2.11 

3.71 

1.60 

2.80 

5.7643286 

1.6510 

MgO 

1.32 


2.39 

2.21 

9.5957763 

0.57470436 

CaO 

1.07 


2.64 

1.99 

(11.4229222) 

0.444441 

SrO 

0.73 


2.98 

1.65 

13.5325857 

0.33798 

BaO 

1.03 


2.68 

1.95 

14.3325535 

0.3126140 

BeS 

2.11 

2.63 

0.52 

2.36 

7.0304618 

0.79059 

MgS 

1.32 


1.31 

1.86 

(13.7042728) 

0.2136* 

CaS 

1.07 


1.56 

1.68 

17.7617688 

0.1766757 

SrS 

0.73 


1.90 

1.39 

(23.444936) 

0.1328* 

BaS 

1.03 


1.60 

1.65 

(25.954702) 

0.10331 

SrSe 

0.73 

2.57 

1.84 

1.37 

(41.860451) 

0.0592* 

IIIA oxides and chalcogenides 






BO 

2.23 

3.71 

1.48 

2.88 

6.52093998 

1.7820 

AIO 

1.48 


2.23 

2.34 

10.0419499 

0.64136 

GaO 

1.77 


1.94 

2.56 

12.982250 

[(0.4271)] 


Contd.. Table 1 
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Contd.. Table 1 


InO 

1.36 


2.35 

2.25 

14.0404513 

0.3296* 

BS 

2.23 

2.63 

0.40 

2.42 

8.1893677 

0.79489 

A1S 

1.48 


1.15 

1.97 

14.6327874 

0.2799 

AlSe 

1.48 

2.57 

1.09 

1.95 

(20.171277 4 ) • 

0.1528* 

IV A oxides and chalcogenides 






CO 

2.73 

3.71 

0.98 

3.18 

6.58620870 

1.93128075 

SiO 

1.87 


1.84 

2.63 

10.1767074 

0.7267521 

GeO 

2.01 


1.70 

2.73 

13.1496247 

0.4856981 

SnO 

1.59 


2.12 

2.43 

14.112331 

0.3557190 

PbO 

1.67 


2.04 

2.49 

14.8526391 

0.30730 

CS 

2.73 

2.63 

0.10 

2.68 

8.72519425 

0.8200434 

SiS 

1.87 


0.76 

2.22 

14.9206887 

0.30352788 

GeS 

2.01 


0.62 

2.30 

22.318826 

0.1865676 

SnS 

1.59 


1.04 

2.04 

25.241415 

0.1368619 

PbS 

1.67 


0.96 

2.10 

27.7119406 

0.11632307 

CSe , 

2.73 

2.57 

0.16 

2.65 

(10.43336091) 

0.5750 

SiSe 

1.87 


0.70 

2.19 

20.7224688 

0.1920117 

GeSe 

2.01 


0.56 

2.27 

(38.401003) 

0.0899* # 

SnSe 

1.59 


0.98 

2.02 

47.954285 

0.06499776 

PbSe 

1.67 


0.90 

2.07 

57.732418 

0.05059953 

SiTe 

1.87 

2.16 

0.29 

2.01 

(23.019425) 

0.16998* 


Contd.. Table 1 
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Contd.. Table 1 


GeTe 

2.01 


0.15 

2.08 

47.11256 

0.06533821 

SnTe 

1.59 


0.57 

1.85 

(62.35 1952 6 ) 

0.4247917 

PbTe 

1.67 


0.49 

1.90 

(79.96 1030 2 ) 

0.03130774 

IVB oxides and chalcogenides 






TiO 

1.42 

3.71 

2.29 

2.30 

11.9938851 

0.53541 

ZrO 

1.29 


2.42 

2.19 

13.5790589 

0.42263 

HfO 

1.58 


2.13 

2.42 

(14.68922344) 

[0.3860] 

TiS 

1.42 

2.63 

1.21 

1.93 

(19.1816180) 

0.1823* 

VA oxides and chalcogenides 






NO 

3.26 

3.71 

0.45 

3.48 

7.46643320 

1.67195 

PO 

2.21 


1.50 

2.86 

10.5479381 

0.7337 

AsO 

2.31 


1.40 

2.93 

13.1809338 

* 0.48482 

SbO 

1.79 


1.92 

2.58 

(14.1261070 

0.3850*** 

BiO 

1.77 


1.94 

2.56 

14.8577 

0.3034 

NS 

3.26 

2.63 

0.63 

2.93 

9.7380290 

0.769602 

PS 

2.21 


0.42 

2.21 

15.7325010 

(0.29) 

BiS 

1.77 


0.86 

2.16 

27.7296877 

[0.112764] 


Note : 


1. Uncertain quantities in the data as mentioned in reference 10 are given in parentheses. 

2. Quantities in square brackets [ ] in column 7 refer to B 0 the value of the rotational constant for 

the lowest vibrational level. 

3. We have considered only the magnitudes of Ax, the electronegativity difference. 

4. The values of B e indicated by *sign are obtained by us by extrapolation from the individual 
graphs. 

5. # Value of B e for GeSe reported as 0.09634051 cm 4 for l X + state in reference 11. 

6. ^Value of B e for SbO reported as 0.3580 cm' 1 for 2 3 4 5 6 FIi /2 state in reference 11. 
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Fig. 1— Plot showing the variation between the values of log B e and 
log (Ax/f* 2 ) for IIA oxides and chalcogenides. 



Fig. 2 -Plot showing the variation between the values of log B e and 
log (Ax/n 2 ) for IVB oxides and chalcogenides. 
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Discussion 

IIA and IVB oxides and chalcogenides 

It can be seen from Fig. 1 that except for BeO and BeS, all other IIA oxides and 
chalcogenides fall on the straight line with m = 0.68 and C = 0.79, and the coefficient 
coefficient R 2 = 0.9864. This suggests that the chemical bonding in these molecules 
(except BeO and BeS) is more or less of the predominantly ionic type, since these 
molecules get systematised with the help of equation (2). The scatter in the case of 
BeO and BeS will be discussed later in this paper. We noticed that the data on 
rotational constant B e of MgS, SrS and SrSe molecules are not available in literature. 
Hence we have used the curve in Fig. 1 to estimate the value of B e in these molecules, 
which come out to be 0.2136 cm' 1 , 0.1328 cm' 1 and 0.0592 cm' 1 respectively. 

In the case of IVB oxides and chalcogenides (Fig. 2), the curve of log B e versus 
log (A^/p 2 ) is found to be a straight line with m = 0.66, C = 0.89 and R 2 = 0.8743, 
passing through the points for TiO, ZrO and HfO. Although the point for TiS lies 
away from the points of the other three oxides, we have ventured to estimate the value 
of B e of TiS from this curve by extrapolation which comes out to be 0.1 823 cm' 1 . 

The other oxides and chalcogenides 

The scatter in the plots of log B e vs. log (Ax/p 2 ) for IIIA, IV A and VA oxides and 
chalcogenides (plots not shown here, as mentioned earlier) indicates that the bonding 
character in these molecules is not predominantly ionic. Hence we thought it 
worthwhile to re-examine the bonding picture in these molecules. The bonding in the 
IV A oxides and chalcogenides may be considered as a very typical case where in the 
molecule AB, the atom A can be thought to be donating two electrons to the atom B 
belonging to VIA group of elements. 

As shown in Fig. 3 (a) the outermost electronic configuration of group IV A 
elements (C, Si, Ge, Sn, Pb) is ns 2 np~ (atom A) and that for group VIA elements (S, 
Se, Te) is ns 2 np A (atom B). In the ionic picture of the molecule (Fig. 3(b))one may 
expect that two p electrons from atom A are transferred to the p orbital of the atom B, 
thereby making the electronic configuration of A 2 ' as ns 2 and that of B 2 as ns 2 np 6 . But 
because of the large scatter in the log B e vs. log (Ax / p 2 ) plot for these molecules, we 
decided to look for an alternative bonding scheme. 
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ns 2 np 6 


/[ A f 3(c) 

I—— ! > Alternate 

ns 2 np 3 possible 

scheme 
suggested 
by the 
ns 2 np 3 authors 




Fig. 3-Bonding schemes in IV A oxides and chalcogenides. 



loo ((XaXb) 1 'V) 












*OQ((ZAXB) m /^) 

Fig. 6-Piot showing the variation between the values of log B e and 
log B / i-i 2 )for VA oxides and chalcogenides. 
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As shown in Fig. 3(c), it is possible that a p electron from the atom B is 
transferred to the p orbital of the atom A, thereby making the electronic configuration 
of A' and B + similar (ns 2 np 3 ), which may lead to electronegativity equalization and in 
turn to the formation of a co valent bond in the molecule. We have checked this idea 
by plotting log B s vs. log ( /p 2 ) for the IV A oxides and chalcogenides. The 
plot thus obtained is shown in Fig. 4 which is a good straight line with m = 0.71 and C 
= 0.94, and correlation coefficient R 2 = 0.9884. This confirms our idea of the bonding 
picture A'B + arising out of electronegativity equalization leading to prominent 
covalent character of the bond instead of the ionic picture A 2+ B 2 ' in IV A oxides and 
chalcogenides. From the plot in Fig. 4 we have estimated the values ofi? e for SiTe and 
GeSe which come out to be 0.16998 cm' 1 and 0.0899 cm' 1 respectively. 

The log j B e vs. log (Ax/p. 2 ) plots for IIIA and VA oxides and chalcogenides also 
exhibit large scatter suggesting that the bonding in these molecules may not be 
predominantly ionic. We therefore decided to check whether the bondin g in the m is 
predominantly covalent. We have plotted the values of log B e vs. log ( fx/ig /p 2 ) 
for IIIA and VA oxides and chalcogenides in Figs. 5 and 6 respectively. As can be 
seen, the plots are good straight lines with negligible scatter. From the plot in Fig. 5 
which is a straight line with m = 0.90, C = 1 .26 and R 2 = 0.983, we have estimated the 
values of B e for InO and AlSe which come out to be 0.3296 cm' 1 and 0.1528 cm' 1 
respectively. The plot for VA oxides and chalcogenides (Fig. 6) is a straight line with 
m = 0.87, C = 1 .23 and R 2 = 0.9893. From this plot we have estimated the value of B a 
for SbO which comes out to be 0.3850 cm' 1 . 

The case of B e O and BeS 

The case of BeO and BeS has been discussed in our recent paper 16 . The scatter for 
the points for BeO and BeS in the plot in Fig. 1 indicates that these molecules may 
have different bonding characteristics as compared to the other oxides and 
chalcogenides of the same group. In the latter molecules, the bonding seems to be 
predominantly ionic as they can be systematised very well with the help of equation 
(2). Hence we decided to check whether these two molecules exhibit covalent 
ch aracter . Shown in Fig. 1 by § are the values of log B e plotted against log 
(s[XaXb / P- 2 ) for BeO and BeS, which are seen to be better accommodated on the 
curve. This confirms our feeling that in BeO and BeS the bonds may have a large 
covalent character. The case of these two molecules is interesting and needs further 
investigation. The predominant covalence in these molecules may probably arise from 
the lightness and small size of the.beryllium atom. 
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Support to the values ofB e estimated by us 

We may mention here that no claim can be made about the accuracy of the values 
of B e estimated by us since the values of the reduced mass in all these molecules 
(except for In O) are reported to be uncertain (see the footnote of Table 1) and also 
because the values of X are known up to two places after the decimal point. 

As is evident from equation (1), the value of the intemuclear distance r e can be 
calculated from the value of B e . To check the correctness of the values of B e estimated 
by us, we have calculated the values of r e in all these nine molecules by using 
equation (1). Given in Table 2 are the values B a estimated by us and those of r e 
calculated from these values of B e . In column 4 of this table are given the 
experimentally determined values of r e in some of these molecules as reported in the 
CRC Handbook of Chemistry and Physics 11 ' 17 . We find that the percentage error of the 
calculated values of r e compared with the experimental values is about five percent. 
Thus, there appears to be a fair agreement between the values of r e determined from 
the B e values estimated by us and the experimental values. Our values of r e can be 
considered to be accurate up to 0.01 A 0 . We may mention here that this much 
accuracy is considered sufficient for most physical and chemical work. 


Table 2-Comparison between the values of intemuclear distance r e calculated by us, and the 
experimentally reported values of r e . 


Molecule 

B e (cm -1 ) estimated 
by us 

Value of r e (in A 0 ) from 
the value of B e 

Experimental value of 
r e (in A 0 ) 

MgS 

0.2136 

2.40 


SrS 

0.1328 

2.33 

2.4405* 

SrSe 

0.0592 

2.61 


InO 

0.3296 

1.91 


AiSe 

0.1528 

2.35 


GeSe 

0.0899 

2.21 

2.1346* 

SiTe 

0.16998 

2.08 


TiS 

0.1823 

2.20 

2.0825* 

SbO 

0.3850 

1.76 

1.826** 


♦values from ultraviolet (UV) spectroscopy 
# values from microwave (MW) spectroscopy 
Value for n ]/2 state as mentioned in reference 11. 
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Concluding Remarks 

Our work has once again brought out the usefulness of electronegativity in the 
systematization of the rotational behaviour of molecules. We have been able to 
estimate the indicative values of rotational constant B e for several molecules for which 
data are not available. Support for our values of B e in three cases comes from 
experimental r e values. It should be worthwhile to determine the values of B e for all 
these molecules experimentally. 

We have been able to suggest a bonding scheme in the case of IV A oxides and 
chalcogenides. More rigorous theoretical calculations for understanding the nature of 
bonding in these and other molecules seem to be necessary. In particular, it will be 
interesting to calculate the exact number of electrons transferred from the atom A to 
the atom B in the diatomic chalcogenides, as has been done for the IA halides by 
Ohwada 18 . 
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